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We review recent progress in attaining a quantitative un- 
derstanding of the scarring phenomenon, the non-random be- 
havior of quantum wavefunctions near unstable periodic or- 
bits of a classically chaotic system. The wavepacket dynam- 
ics framework leads to predictions about statistical long-time 
and stationary properties of quantum systems with chaotic 
classical analogues. Many long-time quantum properties can 
be quantitatively understood using only short-time classical 
dynamics information; these include wavefunction intensity 
distributions, intensity correlations in phase space and cor- 
relations between wavefunctions, and distributions of decay 
rates and conductance peaks in weakly open systems. Strong 
deviations from random matrix theory are predicted and ob- 
served in the presence of short unstable periodic orbits. 



I. INTRODUCTION 

Wavefunction scarring is the anomalous enhancement 
of quantum eigenstate intensities along unstable periodic 
orbits of a classically chaotic system; it is surely one of 
the most visually striking properties of quantum chaotic 
behavior. Observed numerically in unpublished work by 
McDonald (l), scars were later brought to the attention 
of the physics community by Heller ||^, who also pro- 
vided the first theoretical explanation for their existence. 
Numerical evidence and associated analytical work (fol- 
lowed later by experimental tests in a variety of systems) 
showed that scarring was a statistically significant correc- 
tion to Berry's early conjecture ||^ that eigenstates of a 
classically ergodic system are evenly distributed over the 
energy hypersurface, in the semiclassical limit ?i — > 0. 

The structure of quantum eigenstates of a system 
whose classical analogue is ergodic has long been a mat- 
ter of interest for physicists and mathematicians alike. 
When an integrable system is quantized, the resulting 
eigenstates can be well understood in terms of the WKB 
quantization condition along the classical periodic or- 
bits For classically ergodic systems, such a direct cor- 
respondence between quantum wavefunctions and classi- 
cal periodic orbits is, however, not possible. Because 
the typical classical trajectory evenly explores all of the 
available energy hypersurface, it is natural to suppose 
that the quantum eigenstates also have constant inten- 
sity over this entire hypersurface, up to the inevitable 
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quantum fluctuations. Rigorous results have been ob- 
tained on one important aspect of this problem. Theo- 
rems by Schnirelman, Zelditch, and Colin de Verdiere j|] 
consider a classically defined operator and show that its 
expectation value over almost all individual eigenstates 
converges (in the semiclassical limit) to the ergodic, mi- 
crocanonical average of the classical version of the oper- 
ator. In this limit, however, taking the expectation value 
means averaging over more and more de Broglie wave- 
lengths of the wavefunction, since the classical symbol of 
the operator in phase space is assumed to be kept fixed 
as ?i — > 0. Thus, these ergodicity theorems provide lit- 
tle information about the structure of eigenstates at the 
single-wavelength scale (or more generally on scales of a 
single momentum-space channel or over cells of size 0{h) 
in phase space). 

Classical ergodicity by itself is certainly not sufficient 
to ensure ergodically distributed quantum wavefunctions 
over individual wavelengths or channels. Thus, for ex- 
ample, one can consider the "slow ergodic systems," |^ 
where the classical exploration of phase space is suffi- 
ciently slow that the corresponding quantum dynamics 
can only explore a fraction of the available quantum 
channels by the Heisenberg time [this being h over the 
mean quantum level spacing, the time at which the quan- 
tum exploration of phase space effectively ends and the 
quantum dynamics becomes quasi-periodic] . In such sys- 
tems, the participation ratio, i.e. the fraction of available 
quantum channels with which a typical eigenstate has 
significant overlap, approaches zero in the semiclassical 
limit ?i — > 0. However, for any given high-energy eigen- 
state, this small fraction of bright channels are evenly 
distributed over the total phase-space, so that eigenstate 
intensity integrated over a classically-defined patch of 
phase space does in fact converge to the classically re- 
quired value, in perfect accordance with the results of 
Schnirelman, Zelditch, and Colin de Verdiere (SZCdV). 

Such anomalies are not known to arise for strongly 
chaotic systems, by which wc mean systems that have 
everywhere positive entropy. In such systems, the mix- 
ing time (the time for a phase space cell of size Planck's 
constant h to explore all of the available phase space on 
a mesh of that scale) scales only logarithmically with h, 
in contrast with the Heisenberg time, which scales as a 
power law. Thus, in the limit of small h a given classi- 
cal system will always have a mixing time much shorter 
than the Heisenberg time, and many trajectories will ex- 
ist linking any /i-sized cell with any other well before in- 
dividual eigenstates begin to be resolved. It is reasonable 
to conjecture that for such systems individual eigen- 
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states will in fact evenly cover all of the available energy 
hypersurface, up to uncorrelated Gaussian fluctuations 
given by random matrix theory (RMT). Indeed, RMT 
does correctly predict much of the structure of eigen- 
states (and also the eigenvalue spectrum) of classically 
chaotic systems. 

However, we also know that the Gutzwiller trace for- 
mula 1^ is quite effective in using short periodic orbits to 
predict non-random fluctuations in the quantum spectra 
of classically chaotic systems. It is then not too sur- 
prising that the unstable periodic orbits of a chaotic 
system should also have a strong effect on the eigen- 
state structure, modifying RMT predictions. In fact, it 
appears that in the h —^ limit, scarring constitutes 
the most important deviation from perfect ergodicity of 
chaotic eigenstates at the ?i-scale, in the sense of RMT. 
Other deviations from RMT, such as dynamical localiza- 
tion, disappear in the semiclassical limit, but the scar- 
ring phenomenon survives. It, along with symmetry ef- 
fects, is sufficient to explain quantitatively many prop- 
erties of quantum chaotic systems, including wavefunc- 
tion intensity statistics, wavefunction and phase space 
intensity correlations, channel-to-channel transport, res- 
onance widths, and conductances. 

We conclude this introduction by defining what we 
mean by scarring and addressing four common miscon- 
ceptions about this phenomenon; our comments here may 
help to clarify issues which in the past have led to con- 
fusion and even controversy surrounding this subject. 

(i) Scarring is defined as the anomalous enhancement 
or suppression of eigenstate intensity on or near an un- 
stable periodic orbit and its invariant manifolds]^, in a 
chaotic system. Of course, stable periodic orbits also at- 
tract eigenstate intensity, as mentioned above, but the 
reasons for this are well understood in terms of the semi- 
classical theory of integrable systems [Q. More impor- 
tantly, the increased tendency of a probability distribu- 
tion launched in the vicinity of a stable orbit to overlap 
with itself at long times can be understood as a purely 
classical phenomenon. Phases are necessary only to ob- 
tain the quantization energies, in accordance with the 
WKB condition. On the other hand scarring can be 
thought of in the time domain as the increased probabil- 
ity of a quantum wavepacket launched near an unstable 
orbit to overlap with itself at very long times. Classi- 
cal dynamics in a chaotic system shows no such behav- 



^The invariant manifolds are the set of phase space points 
that classically approach the periodic orbit as t +00 (the 
stable manifold), and the set that approach the same orbit 
as t ^ —00 (the unstable manifold). In a two-dimensional 
phase space, these are both curves, and the lines tangent to 
the curves at the periodic point are known as the stable and 
unstable directions, respectively. The (infinitely many) points 
of intersection of the stable and unstable manifolds are the 
homoclinic points of the orbit (see below). 



ior: classically an evolved probability distribution loses 
at long times all memory of its initial position. So scar- 
ring is at first sight a paradoxical phenomenon because 
its presence implies that quantum mechanics retains a 
much better memory at long times of short-time classi- 
cal behavior than does long-time classical mechanics it- 
self. This paradox, as we will see, is resolved by realizing 
that scarring is inherently an interference effect, and can 
only be understood by including phase information. 

(ii) Scarring is a statistically significant deviation from 
single-channel quantum ergodicity, as defined by random 
matrix theory (which predicts uniform wavefunction am- 
plitude over all of available phase space, up to Gaussian 
random fluctuations). Of course random enhancements 
in wavefunction intensity occur even in the absence of 
any periodic orbit; they can be observed, for example, 
in random superpositions of plane waves of constant en- 
ergy Thus, anecdotal evidence, in the form of pictures 
of scarred states, is generally not a good way to test a 
theory of scarring, because even in RMT there is a flnite 
probability to observe any given enhancement in one or a 
few states. An exception to this general rule arises when 
the instability exponent A of a given periodic orbit is very 
small. Then as we shall see one expects to find many 
states (more precisely, a fraction of 0(A) of all states) 
having enhancement factors of order 1/A at this orbit, 
and each one of these individually has an exponentially 
small probability to occur within the context of RMT. 
For moderate A, a few states will not produce the desired 
evidence of scarring, and one needs to consider the en- 
tire distribution of eigenstate intensities on the periodic 
orbit, sweeping through energy or some other parameter. 
This distribution can then be compared quantitatively 
with the scar theory prediction (which differs very signif- 
icantly from RMT). 

(iii) Scarring in not the same as the semiclassical eval- 
uation of individual eigenstate intensities; rather it is a 
statistical property of states near a periodic orbit. It is 
true that in certain systems semiclassical methods work 
well all the way to the Heisenberg time, and can be used 
to resolve individual eigenstates [0. However, this re- 
quires putting in at least as much information into the 
semiclassical calculation as there is in the full quantum 
system itself. Exact wavefunctions of individual states 
are sensitive to small changes in system parameters, and 
they provide no intuition about the overall behavior of 
the system. On the other hand, the statistical predic- 
tions of scar theory (concerning the distribution of wave- 
function intensities, mean conductances and conductance 
fluctuations, mean resonance lifetimes, etc.) depend on 
only a few parameters (e.g., the exponent A), and provide 
intuitively useful information about entire ensembles of 
chaotic systems which could not be obtained by using a 
computer to diagonalize a given Hamiltonian. The dis- 
tinction is also very important experimentally because 
the scarring phenomenon is robust under small pertur- 
bations to the system; thus scar theory predictions will 
often be valid in situations where the exact dynamics is 
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not known well enough to allow individual eigenstates 
to be theoretically determined (either quantum mechan- 
ically or semiclassically) . 

(iv) Finally, as explained above, scars are not a threat 
to SZCdV ergodicity; scar theory addresses quantum 
structure at the ?i-level and is to be properly compared 
with RMT, a much stronger condition than the quan- 
tum ergodicity of Schnirelman, Zelditch, and Colin de 
Verdiere. 

The remainder of this paper is organized as follows: 
in the next section we briefly review some relevant ex- 
periments and discuss several theoretical approaches that 
have been used to understand scarring. Because the liter- 
ature is extensive, it will not be possible to do full justice 
to the contributions of the various authors. The bulk of 
the present paper is a survey of results recently obtained 
using the wavepacket dynamics method, a project initi- 
ated by Heller in the seminal paper |^]. Of course, many 
or all of the results presented here can be understood also 
in the context of the other approaches, particularly in the 
context of Gutzwiller periodic orbit theory. However, as 
a unifying framework we choose to use the time-domain 
methodology in which the results were first obtained, and 
make contact with other points of view where appropri- 
ate. 
Thus, 
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HI introduces and reviews the linear 
(short-time) and nonlinear (long-time) theory of scarring, 
as discussed in greater detail in jl^ . The emphasis here 
is on developing the conceptual framework and describ- 
ing the relations among the concepts of the Gaussian 
wavepacket, the short-time and long-time autocorrela- 
tion function, the local density of states (LDOS), the 
linearized energy envelope, and the inverse participation 
ratio (IPR). This machinery is put into use immediately 
in Section where we perform numerical tests and com- 
pare the data with quantitative theoretical predictions 
about spectra, IPR's, and energy correlations. An analy- 
sis of wavefunction intensity statistics follows in Sec- 
tion 0, where we find in accordance with scar theory a 
power-law tail in the intensity distribution, in sharp con- 
trast with the exponential-tail prediction of RMT. Sec- 
tion VI addresses and resolves ambiguities in the measure 



of scarring intensity, leading to a universal and fully opti- 
mal measure in the strong scarring limit p^ ] . Section |VII| 
deals with more recent applications of scar theory to open 
systems, discussing resonances, the probability to remain 
in an open chaotic system , and the effect of scars on 
conductances [^. Finally, in Section VIII wc sum up 
the results and suggest some directions for the future. 



II. EXPERIMENTAL EVIDENCE AND 
THEORETICAL PROGRESS 



A. Experimental and numerical evidence for scars 

Scarred eigenfunctions were first directly observed 
experimentally in two-dimensional microwave cavities 
by Sridhar [|l6|, and independently by Stein and 
Stockman Here one makes use of the exact corre- 

spondence between the quantum wavefunctions of a bil- 
liard system and the electromagnetic modes in a cavity 
with the same boundary conditions. The electromag- 
netic modes at various energies can be mapped either by 
measuring the reflected microwave power as a function of 
energy and position jl^ , or by using a perturbative tech- 
nique involving measurements of small energy shifts [p^ . 
Individual eigenstates scarred by speciflc short unstable 
orbits can be observed, although no attempt was made 
to quantify the scarring strengths. Of course, experi- 
mental limitations may preclude the collection of data 
at extremely high frequencies, where one could test scar 
theory predictions about the statistics of chaotic wave- 
functions in the semiclassical limit. 

Another very important set of experiments |l8| have 
been performed in semiconductor heterostructures. A 
resonant tunneling device contains a quantum well sep- 
arated by a barrier from a two-dimensional electron gas. 
Electrons can be injected into the quantum well, and 
information about electron wavefunctions obtained by 
studying the tunneling rates. When a large tilted mag- 
netic field is turned on, the classical dynamics of the sys- 
tem is known to be chaotic. Indeed, very strongly scarred 
wavefunctions are observed in such systems, associated 
with orbits that have a small instability exponent for a 
wide range of system parameters jiot - Again, the scarring 
patterns are in good qualitative agreement with theoret- 
ical expectations. 

Given the difficulty in collecting large amounts of good 
experimental data at very high energies, numerical sim- 
ulations of quantum chaotic systems have from the very 
beginning been useful in providing "empirical" support 
for the scarring phenomenon. Many different quan- 
tum chaotic systems have been thus explored, including, 
for example, the hydrogen atom in a uniform magnetic 
field 1^ . In recent work by Li and Hu , a plane wave 
decomposition method is used to find wavefunctions near 
the one-millionth state of the stadium billiard, showing 
that strong scarring is still present at such energies. Cer- 
tain orbits, including the diamond orbit, are observed to 
be scarred in a corridor of width scaling as the wave- 
length (rather than as the geometric mean of the wave- 
length and the system size; see below). The reason for 
this is not well understood by the present author, though 
it may be associated with the very special properties of 
the stadium, namely, the presence of marginally stable 
orbits and intermittency. 

A statistical analysis of intensities in highly-excited 
wavefunctions of the stadium and other hard chaotic 
systems would prove very useful; such a detailed anal- 
ysis of medium-energy wavefunctions in phase-space co- 
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ordinates on the billiard boundary is now in preparation 
p2| . A somewhat different point of view is taken by 
Backer, Schubert, and Stifter ||2^, who examine the ap- 
proach to SZCdV quantum ergodicity in position space 
for various chaotic billiards (again, in a medium energy 
range). Simonotti, Vcrgini, and Saraceno have stud- 
ied scar intensity on the boundary of a stadium billiard, 
using a phase-space representation and also considering 
improved test states which are superpositions of several 
phase-space Gaussians (see Section VI). 



B. Several theoretical approaches 

The first theoretical explanation of the scarring phe- 
nomenon is found already in where a Husimi space 
measure of scar strength is used. In this approach (which 
we discuss at much greater length beginning with the 
following section), wavefunction intensity in phase space 
is defined by the squared overlaps of the wavefunction 
with Gaussian wavepackets \a''°'P°) centered at all possi- 
ble phase space points (goiPo)- If we consider one such 
Gaussian and look at its squared overlaps with all eigen- 
states of the system, we obtain a local density of states 
at the location of that wavepacket 
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What Heller realized was that the autocorrelation func- 
tion of the wavepacket. 



(2) 



the Fourier transform of the spectrum S{E), shows large 
recurrences at short times which are integer multiples of 
the orbit period Tp, if the point {qo,Po) is on or near 
a periodic orbit of the system. The strength of these 
recurrences falls off exponentially with time, roughly as 
g-At/Tj=^ where A is the instability exponent of the un- 
stable periodic orbit. In the semiclassical limit ?i ^ 0, 
these periodic recurrences (both amplitudes and phases) 
can be easily analytically computed in terms of the clas- 
sical stability matrix and the classical action of the pe- 
riodic orbit. The recurrences are called "linear", as they 
are obtained by linearizing the dynamics around the un- 
stable periodic orbit in question. This short-time struc- 
ture in A(t) produces an envelope in the spectrum S{E) 
which must be the smoothed version of the true LDOS at 
the periodic orbit. Thus, by considering only the short- 
time dynamics of the Gaussian wavepacket, we obtain an 
energy-smoothed LDOS, 
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from which we can easily extract the average wavefunc- 
tion intensity near various energies. If the exponent A is 
small, the smoothed LDOS is shown to have a sequence 



of bumps separated by 2-kTi/Tp and of width scaling as 
Xh/Tp. At these optimal scarring energies (which are 
analogous to the EBK energies for an integrable sys- 
tem), wavefunctions must, on average, have larger than 
normal intensities on the periodic orbit, while far from 
these spectral peaks wavefunctions must be typically an- 
tiscarred, i.e. have less than expected intensity on the 
same orbit. 

The linear theory of scarring provides information only 
about the average wavefunction intensity in some energy 
range; it does not tell us whether most of that inten- 
sity is contained in a few wavefunctions (the "totalitar- 
ian" scenario) or whether it is equally distributed among 
all the wavefunctions near that energy (the "egalitarian" 
scenario). For this, we need a nonlinear theory of scar- 
ring, one that takes into account long-time recurrences 
in the autocorrelation hmction A{t). These recurrences 
correspond to pieces of the wavepacket leaving the vicin- 
ity of the periodic orbit, traveling through other regions 
of phase space where they are subjected to complicated 
chaotic dynamics, and eventually returning to the peri- 
odic orbit. In many situations, these long-time recur- 
rences may be treated statistically, producing a sepa- 
ration of scales between the classically-obtained short- 
time behavior near the periodic orbit and the random 
long-time behavior, as suggested already by Antonsen et 
al. pa. What is important to note here is that the short- 
time behavior must necessarily leave an imprint on the 
long-time behavior and on stationary properties like the 
LDOS, as we will see more clearly in the following sec- 
tion. 

First, we briefly mention some alternative theoretical 
approaches to the scarring problem, which provide useful 
and complementary perspectives on the subject. A the- 
ory of wavefunction scarring in position space has been 
developed by Bogomolny clearly the position space 
basis is a very relevant one for many experimental ap- 
plications. Bogomolny again smoothes the wavefunction 
intensity over some small energy range AE', 



(4) 



and (in the semiclassical limit ?7, — > 0) represents this 
smoothed intensity as a sum over periodic trajectories of 
the system: 



< 



Im( Bp{z) exp 
p \ 



2h 



VkVn 



(5) 



AE 



For each periodic orbit p, z is chosen to be a coordi- 
nate along the orbit, while the j/m ("t- = 2 . . . n) are all 
the coordinates perpendicular to the orbit. is the 
classical action around the orbit, while is a focusing 
factor which can be easily obtained from the stability 
matrix of the orbit (i.e. from the linearized dynamics 
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for small y,y). Similarly W^™ measures the quadratic 
change in the action for small values of y, as one moves 
away from the orbit. It is important to note here that 
the sum in Eq. ^ is effectively finite: any orbits of period 
Tp >> fi/lS.E can be dropped as their oscillatory con- 
tribution vanishes when integrated over the smoothing 
interval Ai?. As IS.E becomes small, more and more or- 
bits need to be included in the sum (the number of orbits 
grows exponentially with period Tp in a chaotic system) . 
The constant po{q) term can be thought of as coming 
from the zero-time dynamics of the system; it is related 
to the mean density of states at energy E. Bogomolny 
computes the oscillatory contributions to the above sum 
from various short periodic orbits in the stadium billiard 
and finds qualitative agreement with numerical data. 

Bogomolny's semiclassical Green's function approach 
is clearly closely related to our wavepacket dynamics 
method, as the semiclassical Green's function can be ob- 
tained from the semiclassical time-domain propagator by 
a stationary-phase Fourier transform. One difference be- 
tween the approaches is that Bogomolny envisions sum- 
ming over a large number of periodic orbits to get as close 
as possible to an energy domain resolution of order of a 
mean level spacing. As mentioned above, in some sys- 
tems it is indeed possible to use semiclassical methods to 
compute individual eigenstates of the system jl^ . In fact 
for this purpose one needs information only about orbits 
of period up to the mixing time (which scales logarith- 
mically with Ti) rather than the Heisenberg time (which 
scales as a power law). However, our aim here is to make 
predictions about the distribution of scarring strengths 
based only on linearized information around one periodic 
orbit; for this purpose most other orbits which produce 
additional oscillations in the density of states may be 
treated statistically |25|. It is important to note in this 
context that if we are measuring wavefunction intensities 
on a given orbit p of period Tp, then in the semiclassical 
limit there are no short orbits that come close to this or- 
bit (on a scale of K) in phase space. The only oscillatory 
contributions which will need to be taken into account 
arc from orbits closely related to orbits homoclinic to p 
[homoclinic orbits are those that approach p at large neg- 
ative times, perform an excursion away from p into other 
regions of phase space, and then again approach p at large 
positive times]. In fact, in the ?i — > limit the periodic 
orbit sum (Eq. ^ for a point g on a given periodic or- 
bit p can be written equivalently as a contribution from 
orbit p itself plus a sum over trajectories homoclinic to 
p. Although the two points of view are mathematically 
equivalent, the homoclinic sum approach makes explicit 
the special role of the orbit p near which we are making 
measurements . 

We also take this opportunity to note that a position 
space basis, though obviously physically natural in many 
measurement situations, is not generally an optimal one 
for detecting scar effects. Unless the periodic point q 
happens also to be a focusing point of classical trajec- 
tories near the orbit, only a small fraction of the total 



scar strength is captured in the position basis, and the 
fraction becomes smaller as h decreases (or as the en- 
ergy increases) . An easy way to see this is to notice that 
the effects of a classical trajectory in quantum mechan- 
ics generically extend to a region around the orbit scaling 
not as a wavelength but rather as the square root of a 
wavelength (and similarly the affected region scales as the 
square root of the total number of channels in momentum 
space). Thus, unless either the stable or unstable mani- 
fold of the orbit p at periodic point q happens to be ori- 
ented along the momentum direction, the position space 
basis will not be optimal, as reflected in the falling off of 
the focusing prefactor Bp with energy (and similarly the 
momentum basis will not be optimal, unless one of the 
two invariant manifolds is oriented along the position di- 
rection) . All this will become more clear in the exposition 
of the following section. In any case, one should keep in 
mind that a position space basis can always be considered 
as a special limiting case of the Gaussian wavepacket test 
state, where the position uncertainty of the wavepacket 
becomes comparable to a wavelength, and the momen- 
tum uncertainty becomes large. 

Wigner phase space analysis of the scarring phe- 
nomenon was pioneered by Berry p7| . Berry consid- 
ers the Wigner spectral function, W{x, E,AE), again 
smoothed over an energy interval AE near E. x = 
{qo,Po) is a phase space point. Being formulated in phase 
space, the approach more closely resembles that of Heller, 
and working in Wigner phase space instead of Husimi 
space also eliminates the need to choose the (apparently 
arbitrary) eccentricity and orientation of the Gaussian 
wavepackets a'^" . The downside of Wigner phase space 
is the absence of a positivity condition on W\ thus the 
value of the spectral function cannot be considered as 
corresponding to an intensity or a probability of being 
found near a certain point x (and random matrix theory 
is therefore not applicable). The Husimi function, which 
is manifestly positive definite, can be thought of as a 
phase space smoothing of the Wigner distribution over a 
phase space region scaling as h. The ambiguity in choos- 
ing the Gaussian centered on x over which this smoothing 
is to be performed is indeed an important issue, to be con- 
sidered carefully towards the end of the following section. 
There we will see that to obtain the maximal scarring ef- 
fect, the Gaussian must be chosen to be properly oriented 
along the stable and unstable directions at the periodic 
point. [An arbitrarily large wavepacket width is allowed 
along either of these directions, with a correspondingly 
small width in the orthogonal direction. Also, strong, 
but non-maximal, scarring will generally be obtained for 
any wavepacket with width scaling as ^/h in both the po- 
sition and momentum directions.! See also Section VI for 



a discussion of even better measures of scarring possible 
if one is willing to go beyond the Gaussian wavepacket 
approximation. 

Berry studies carefully the structure of the spectral 
Wigner function as a; moves off the closed orbit within 
the energy hypersurface and also as x moves off the en- 
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ergy surface itself (the fringes he obtains disappear upon 
Husimi averaging). In particular, as x moves off the pe- 
riodic orbit, the transverse oscillations take the form of 
a complex Gaussian, with fringe spacing scaling as Vh. 
The amplitude of the fringes depends only on the stabil- 
ity matrix, although the precise pattern changes as one 
moves along the orbit. In a two-dimensional phase space, 
the fringes form hyperbolas approaching the stable and 
unstable manifolds of the orbit. 

A common limitation of the analyses |H,|2^,^ is that 
they make no predictions about the properties of the 
spectral fluctuations on scales much smaller than h/Tu, 
where Td ^ Tp/X is the decay time of the unstable or- 
bit being studied. Therefore is it not possible to make 
quantitative predictions about individual wavefunction 
intensities, participation ratios, etc., without explicitly 
doing a Gutzwiller sum over all periodic orbits. Even if 
the sum can be performed, it is by no means clear that it 
will converge in all cases (e.g. in systems where caustics 
are important [^). When the sum does converge it may 
produce individual semiclassical wavefunctions very dif- 
ferent from the quantum eigenstates, due to diffraction 
and other "hard quantum" effects. Furthermore, as we 
discussed in the introduction, such Heisenberg-time cal- 
culations are extremely sensitive to small perturbations 
on the system. What one would like is to be able to 
say precisely how often a given singie-wavefunction scar 
strength will appear on a given orbit, at what energy, 
and at what parameter values. In the semiclassical limit, 
this can in fact be done using only information about lin- 
earized dynamics near the orbit itself, and, in some cases, 
about a few strong isolated homoclinic recurrences which 
cannot be treated statistically. 

Before proceeding to the main part of the paper, we 
mention some of the important contributions to scar the- 
ory since the early work of Heller, Bogomolny, and Berry. 
Many different aspects of the scarring phenomenon have 
been theoretically investigated, and we here mention only 
a few which are most relevant to the present work. Agam 
and Fishman |2|] define the weight of a scar by integrat- 
ing the Wigner function over a narrow tube in phase 
space, of cross-section U, surrounding the periodic or- 
bit. The Fredholm method can be used to obtain a semi- 
classical formula for scars. Alternatively, de Polavieja, 
Borondo, and Benito construct a test state highly 
localized on a given periodic orbit using short-time quan- 
tum dynamics. We defer all discussion of improved mea- 
sures of scarring beyond the simple Gaussian wavepacket 
measure until Section where this important issue is 
treated in full generality. 

Klakow and Smilansky have used a scattering ap- 
proach to quantization to study the wavefunctions of bil- 
liard systems. They treat carefully the wavefunctions on 
the Poincare surface of section, and relate their properties 
to scarring in configuration space, de Almeida |Q uses 
the Weyl representation to establish connections between 
classical and quantum dynamics, with particular appli- 
cation to the semiclassical Wigner function and scars. 



Tomsovic pd^ has used parametric variation as a new 
method for studying scar effects; scars are shown to in- 
duce correlations between wavefunction intensities on a 
periodic orbit and the level velocities of these wavefunc- 
tions when certain system parameters are varied. We also 
mention the work of Arranz, Borondo, and Benito ]3^ ] 
who have probed the intermediate region between regular 
and strongly chaotic quantum behavior, and have shown 
how scarred states first arise from the mixing of pairs of 
regular wavefunctions as h is decreased (but well before 
one reaches the semiclassical limit which is the main fo- 
cus of the present review) . Finally, several groups |35|-|37|] 
have studied the hyperbolic scar structures associated 
not only with the periodic orbit itself but with its invari- 
ant manifolds and homoclinic orbits. This subject we 



also return to in Section VI 



III. LINEAR AND NONLINEAR THEORY OF 
SCARRING 

Consider an arbitrary (unstable) periodic orbit of a 
chaotic system. We will simplify the initial exposition by 
taking the orbit to be a fixed point of a discrete-time area- 
preserving map on a two-dimensional phase space. [The 
case of a higher-period orbit of such a map, or of an orbit 
of a continuous-time dynamics in two spatial dimensions, 
can be reduced to the present case simply by iterating 
the original map, or by taking a surface of section map, 
respectively. The generalization to these two situations 
is straightforward, and is discussed towards the end of 
this section and in more detail in [ pl| . Similarly, the 
entire framework generalizes with minor modifications to 
unstable orbits in higher-dimensional systems.] 

The fixed point is taken to be at the origin of phase 
space, {q,p) = (0,0). Without loss of generahty, we can 
take the stable and unstable directions at the fixed point 
to be vertical (p) and horizontal (q), respectively. The lo- 
cal dynamics around the periodic orbit can always be put 
into this form via a canonical transformation of the co- 
ordinates. Then the only parameter describing the local 
(linearized) dynamics near the orbit is A, the instabil- 
ity exponent for one iteration of the orbit. Locally, the 
equations of motion are given by 



P 



q 

P 



e''q 
e-^'p. 



(6) 



We are interested in studying the behavior of the quan- 
tum wavefunctions near this fixed point; as discussed ear- 
lier the natural way to do this is by examining their over- 
laps with a phase-space Gaussian centered on this point: 



aa{q) = 



1 



1/4 



(7) 



This is a minimum-uncertainty state centered at the ori- 
gin of phase space, with width aVh in the g— direction 
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and y/h/a in the p— direction. The width a is an arbi- 
trary parameter: is the aspect ratio of the phase-space 
Gaussian, typically chosen to be of order unity. Any real 
value of cr, as we shall soon see, produces a wavepacket 
capable of optimally measuring scar strength. The ap- 
parent ambiguity in the choice of a is an important issue 
that we will return to in Section [v| . 

In fact, there is another, less obvious, ambiguity in the 
choice of the test wavepacket in Eq. |^. A non-circular 
Gaussian can be made to have any desired orientation 
in the q — p plane by complexifying cr; this produces a 
test state in which q and p are correlated. Any Gaus- 
sian which is not oriented along the stable and unstable 
axes {p and q, in our example) will be non-optimal in 
the sense of having smaller short-time recurrences in its 
autocorrelation function A{t) compared with an optimal 
wavepacket, and therefore having also smaller fluctua- 
tions in its spectral intensities |(n|a)p. A wavepacket 
defined to be circular in a coordinate system in which 
the two invariant directions are not orthogonal will be 
transformed into such a "tilted" wavepacket in our pre- 
ferred coordinate system (Eq. |^) defined by the stable 
and unstable directions. Also, a position or momentum 
state will generically appear "tilted" in this preferred co- 
ordinate system of the periodic orbit. See Eq. ^ below, 
and the associated discussion. 

For small enough h, the wavepacket \a^) and its short- 
time iterates are contained well within the linear regime. 
As long as the wavepacket stays in the phase space region 
surrounding the periodic orbit in which the linearized 
equations of motion (Eq. |6|) apply, the evolution of the 
wavepacket is completely semiclassical, given simply by 
the stretching of the q— width parameter a (and the as- 
sociated shrinking of the momentum- width (Tp = f ^^)- 
More explicitly, at short times we have 



uL 



(8) 



where U is the unitary operator implementing the quan- 
tum discrete-time dynamics, Unn represents the quanti- 
zation of the linearized behavior near the periodic orbit, 
and t is time, measured in units of a single mapping. 
Here —0 is a phase associated with one iteration of the 
periodic orbit: it is given by the classical action in units 
of h, plus Maslov indices if appropriate. 

Then the short-time autocorrelation function of the 
wavepacket is easily seen from Eqs. to have the form 



-i(pt 



(9) 



y/cOSh{Xt) 

The 'lin' subscript indicates that Eq. || describes the 
piece of the autocorrelation function coming from the lin- 
earized dynamics around the periodic orbit. For a weakly 
unstable orbit (small A), ^iin(i) is slowly decaying, with 
strong recurrences happening for the first 0(1/A) itera- 
tions of the orbit. We note that the short-time autocor- 
relation function j4ii„(t) is c-independent (for real cr), a 
fact that will prove important later on. 



After a certain time, the wavepacket leaves the lineariz- 
able region and nonlinear recurrences begin to dominate 
the return probability. This time scale, called the log- 
time, scales logarithmically in N ^ because of the 
exponential divergence of trajectories away from the or- 
bit: 



A 



(10) 



Here N is the total number of /i-sized cells in the ac- 
cessible phase space (also equal to the dimension of the 
effective quantum mechanical Hilbert space), and / is the 
fraction of this phase space (typically 0(1)) in which the 
linearized equations of motion (Eq. |6|) apply. Semiclassi- 
cally, the long-time recurrences are given by a sum over 
all trajectories homoclinic to the original periodic orbit. 
These are classical trajectories that ast—f ~oo approach 
the fixed point along the unstable manifold (the g-axis, 
in our case), and as t ^ -|-oo approach it again, this 
time along the stable manifold (p) . Let us choose a point 
X — {Sq, 0) on the unstable axis which is well within the 
linearizable region surrounding (0,0), and which after T 
time steps maps to a point /^(a;) — {0,Sp) on the sta- 
ble axis, again inside the linearizable regime. In between 
the trajectory must undergo some complicated nonlinear 
behavior in phase-space regions far from the fixed point. 

The path of the full trajectory and its surroundings 
naturally separates into three stages: (i) ti steps during 
which a thin vertical strip centered on x of the original 
wavepacket shrinks vertically and stretches horizontally 
as its center moves out horizontally at an exponential rate 
along the unstable manifold, finally reaching the edge^ of 
the linearizable region; followed (ii) by complicated non- 
linear dynamics for T — ti — T3 steps, which eventually 
brings the center of the distribution back into the lin- 
earizable region; and finally (iii) a time of T3 steps during 
which the dynamics takes us back through the lineariz- 
able region, this time along the stable direction. After T 
steps, the original thin vertical strip has become a long 
horizontal strip which then intersects the original Gaus- 
sian. The total contribution to the wavepacket autocor- 
relation function at time T coming from this homoclinic 
excursion is given by a product of five factors: 



• Q(T-Ti-r3)e^'^" 



(11) 



The factors e ^'^^/'^ and e -^^3/2 are instability factors 
associated with the linearized motion of the wavepacket. 



■^Of course, the exact boundaries of the linearizable region 
are arbitrary and in no way affect the final result. All that 
is required is for the wavepacket itself to be well-contained in 
the region where the linearized equations of motion hold. 
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while e ^'^^'^ and e ^'^^'^ are the corresponding phases. 
The suppression factors e-'^^'/^fia^ ^^^^ ^-Sp^/2hal 
suit from the fact that the initial and final points of 
the excursion are both off-center relative to the Gaus- 
sian wavepacket (recall ap = 1/cr). The total correlation 
function A{T) is given semiclassically by a sum of terms 
of the form (Eq. over all the homoclinic excursions of 
length T: 

v4sc(r) = ^A^,c. (12) 

HC 

Because the long-time homoclinic orbits come back 
with complicated accumulated phases ^noniin, and the 
number of these recurrences grows exponentially with 
time, one might expect the total long-time return am- 
plitudes to be given by Gaussian random numbers. In 
fact, however, we must consider together contributions 
from all homoclinic excursions which lie on a single ho- 
moclinic orbit, e.g. the trajectory which takes x — *■ {x) 
in T steps and the one which takes /^^ {x) f'^^^'^ {x) 
in T + A2 — Ai steps, where Ai and A2 are small 
enough that the resulting phase-space points are still in 
the linearizable regime. These paths all come back in 
phase with each other, having the same nonlinear phase 
'/'noniinj and givc risc to short-time correlations in A(t) 
for large t [Q. These correlations are shown to be re- 
lated to the short-time dynamics of the original Gaussian 
wavepacket. In fact, we can write the return amplitude 
at times Tiog < T < Tf/ {Tr = N is the Heisenberg 
time, where individual eigenstates begin to be resolved) 
as a convolution 

A{T) = J2 A„d(T) An(T - r) . (13) 

T 

Here A^^ is the short-time return amplitude, and A^^d 
has the statistical properties of an uncorrelated random 
Gaussian variable. In particular, 

< Arnd{T) > = 
< A:Ut)A,.^^{t') > =j^Srr'. (14) 

The prefactor 1 /N provides the proper classical normal- 
ization: in the absence of interference effects, the prob- 
ability to come back is equal to the probability for vis- 
iting any other state in the Hilbert space. The average 
in Eq. ^ is taken over long times r, Tiog <C r <C Ty, 
and/or over an ensemble of systems which all have the 
same linearized dynamics around our chosen periodic or- 
bit. In either case, the total size of the Hilbert space N 
(— 1/h for a phase space area normalized to unity) has 
been assumed to be large. We then obtain 

< A{t) > 
< A*{t)A{t + A)> 1 ^ AL(s)An„(s + A) . (15) 



At times beyond the Heisenberg time, this gets modified 
to 

< A*{t)A{t + A) ^ ^ AUs)Aun{s + A) . (16) 

s 

F is a factor associated with the discreteness of the eigen- 
states: it is 3 for real eigenstate-test state overlaps and 
2 for complex overlaps. 

The long-time autocorrelation function is thus self- 
correlated on a scale A ^ X^^ . Qualitatively, this can be 
understood on a purely classical level: once probability 
happens to come back to the vicinity of a weakly unsta- 
ble periodic orbit, it tends to stay around before leaving 
again. On the other hand, the overall enhancement in 
the total return probability at long times: 

obtained by combining the general expression Eq. |l^ 
with the short-time overlap dynamics of the Gaussian 
wavepacket (Eq. ^ , is fundamentally an interference phe- 
nomenon, and signals a kind of quantum localization. In 
terms of homoclinic orbits, the enhancement arises be- 
cause of the coherent addition of paths of the same length 
T lying on the same homoclinic orbit (setting Ai — A2 
in the discussion following Eq. |2|). 

We now define S{E) to be the Fourier transform of the 
autocorrelation function, 

-, +00 

i— — 00 

For a non-degenerate spectrum, it is easy to see (by in- 
serting a complete set of eigenstates) that this produces 
the line spectrum of Eq. ^ Cutting off the sum in Eq. |l^ 
at iTiog, or equivalently by including only linearized dy- 
namics around the periodic orbit, we obtain the smoothed 
local density of states: 

5iin(i?) = ^An„(t)e^^*, (19) 

t 

an envelope centered at quasienergy E = cj) (see Eq. |^), 
of width 6E ~ A, and of height ~ (a factor of 2tt 
has been inserted into the definition of Sn^ for future 
convenience). E — (p is the analogue of the EBK quanti- 
zation condition for integrable systems; here, because of 
the instability of the orbit, scarred states can live in an 
energy range of 0(A) around the optimal energy. States 
with energy more than 0(A| log A|) away from resonance 
tend to be antiscarred, i.e. they have less than expected 
intensity at the periodic orbit [see Eq. ^ and Fig. ^ . 

Now long-time (nonlinear) recurrences as in Eq. |l^ 
lead to fiuctuations under the short-time envelope in the 
full spectrum S{E). Because these recurrences involve a 
random variable convoluted with the short time dynam- 
ics, in the energy domain we obtain random fluctuations 
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multiplying the short-time envelope. (It is easy to see 
physicaUy that the random oscillations must multiply the 
smooth envelope: if they were merely added to it, the to- 
tal spectrum would go negative away from the peak of 
the envelope.) Finally, at the Heisenberg time Th — N, 
individual states are resolved [0,^, and we see a line 
spectrum (Eq. with a height distribution given by 



Ina^ = |(n|ao-)P = ranSvm{En) ■ 



(20) 



Here ran are random variables (with mean < ran >— 
1 /N) drawn from a chi-squared distribution of one degree 
of freedom (two degrees of freedom for complex {n\aa-))- 
Thus we obtain a random (Porter-Thomas) line spec- 
trum, multiplying the original linear envelope [see Fig. |l| 
for a numerical example] . 

We now need to discuss the concept of an inverse par- 
ticipation ratio (IPR), a very useful measure for studying 
deviations from quantum ergodicity. We define 



(21) 



(Note that J^n = 1 by normalization.) Being the 
first non-trivial moment of the eigenstate intensity (/na„ ) 
distribution, the IPR provides a convenient one-number 
measure of the strength of scarring (or any other kind of 
deviation from quantum ergodicity). The IPR would be 
unity for a wavepacket that had equal overlaps with all 
the eigenstates of the system; the maximum value of TV 
is reached in the opposite (completely localized regime), 
when the wavepacket is itself a single eigenstate. Random 
matrix theory predicts an IPR of F, the strong quantum 
ergodicity factor defined above in Eq. [l^ . 
However, from Eqs. plJlS we see that 



IPR,„ = hm ^5]|A(t)| 



T-l 



T 



(22) 



as one might expect, localization is associated with an 
enhanced return probability at long times. Now from 
Eq. |l^ we see that scar theory predicts an enhancement 
in the IPR over random matrix theory: 



IPR„^ -^^^ 



cosh(As) 



F 



(23) 



where the limit of small A has been taken. {F, as before, 
is 3 or 2, depending on whether the states are real or 
complex, respectively.) The IPR thus decomposes into 
a product of two contributions: the shape of the short- 
time envelope coming from linear dynamics around the 
periodic orbit, and a quantum fluctuation factor F, as 
predicted by Porter-Thomas statistics. See Fig. || in the 
following section for numerical evidence in support of this 
prediction. 

Before concluding this section, we briefly mention sev- 
eral straightforward generalizations of the above analy- 
sis. Eq. ^ for the short-time autocorrelation function A 



and the linear spectral envelope 5'iin which follows from 
it apply to a fixed point of a discrete-time map, but the 
analysis extends easily to higher-period orbits and to dy- 
namics in continuous time. For a period Tp orbit of a 
map, the linearized autocorrelation function of an opti- 
mally oriented wavepacket takes the form 



^lin(t) 



-i4>t/Tp 



y/cosh{Xt/Tp) 



JtmodTp,0 : 



(24) 



where A is now of course the instability exponent for 
one iteration of the entire orbit, and </> the correspond- 
ing phase. This result is the same on each of the peri- 
odic points of the orbit; however, the optimal Gaussian 
wavepacket will in general have a different orientation at 
each point, due to the rotation of the invariant manifolds 
as one moves along the orbit. The resulting smoothed 
local density of states Sun then has Tp peaks in the 
quasienergy interval [0, 27r), each with a height scaling as 
A~^ as before (in the regime of small instability A ^ 1), 
and a width scaling as X/Tp. Notice that the time and 
quasienergy coordinates have been simply rescaled by Tp; 
the ratio of scar peak separation to peak width, as well 
as the inverse participation ratio, remain independent of 
the period Tp, both scaling simply as A~^. 

For a continuous-time Hamiltonian dynamics, the 
spectral envelope has infinitely many peaks, with en- 
ergy spacing h/Tp, height scaling as A~^ as before, and 
width scaling as hX/Tp for small A. The positions of 
the peaks correspond of course to the EBK quantiza- 
tion condition along a stable orbit, and the widths can 
be thought of as resonance widths due to the instability 
of the orbit. In the semiclassical limit, more and more 
states are found under each of the scarring peaks, and the 
individual line heights are distributed in accordance with 



Eq. |2C| (see 
biUiard). T: 



| 22| for a numerical example in the stadium 
"he continuous-time dynamics also contains 
an additional time scale not found in the discrete-time 
system, namely the time Tfroc for the test wavepacket 
to traverse itself near t = and again every time it re- 
turns to the vicinity of its original location. Thus, the 
(5-function sum in Eq. |2j must, for a continuous-time dy- 
namics, be convoluted with a very short-time window of 
width Tfrcc (Tfrco <C Tp in the semiclassical limit). In the 
energy domain, this results in the infinite sequence of 
scarring peaks being multiplied by a very broad envelope 
of width Efj-cc — h/Tfycc- Efrcc is the energy uncertainty 
of the test wavepacket (Sfrcc ^ pVha-p = p\/h/az), 
and evidently depends on , the position-uncertainty of 
the wavepacket in units of ^/K along the direction of the 
orbit. This broad energy envelope is a classical effect (i.e. 
a classical phase-space distribution corresponding to the 
quantum wavepacket has this distribution of energies), 
and must first be taken into account before deviations 
from quantum ergodicity, such as scarring, can be ana- 
lyzed. 



See 
these issues. 
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for a more complete discussion of 



lin 



The individual peak height distribution (Eq. GOl) in the 
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local density of states applies in the semiclassical limit for 
a hard chaotic system. In this limit, the homoclinic (non- 
linear) recurrences are all weak, and no individual recur- 
rence can affect significantly either a particular eigenstatc 
or statistical wavefunction properties such as the IPR. 
However, at finite h it is quite possible for a certain class 
of identifiable homoclinic recurrences at medium times 
Tmed to have a significant effect on the spectrum. The 
contributions from these special recurrences can then be 
computed analytically, and an improved spectral enve- 
lope created which takes into account not only the peri- 
odic orbit dynamics on scale Tp/A, but also the special 
strong recurrences at the time scale Tmcd- The general 
procedure, discussed in more detail in ||ll|, is a sequen- 
tial analysis of the various time scales, beginning with the 
shortest times. Thus, we start with the classical envelope 
coming from the free dynamics on time scale Tftoc, then 
we include the effects of the short periodic orbit on which 
our wavepacket is located, then the effect of any special 
class of strong homoclinic recurrences (an ?i-dependent 
contribution), and finally the effects of random recur- 
rences at long times, all the way to the Heisenberg time 
Th- In [ pT| an example is given of a system where the 
phase-space mixing away from the periodic orbit being 
studied is slow, so that at finite U very distinct oscilla- 
tions in the LDOS at the orbit are obtained which can 
be assigned to a special class of intermediate-time recur- 
rences. 

The expression (Eq. |^) is valid, as we have discussed, 
for any wavepacket optimally oriented along the invari- 
ant manifolds of the periodic orbit. There is a one- 
parameter family of such optimal wavepackets at each 
periodic point: in the coordinate system of Eq. ^jin which 
the stable and unstable directions are orthogonal, this 
free parameter is the width of the Gaussian along either 
of the two directions. Gaussians that are tilted relative 
the invariant directions are non-optimal: the degree of 
non-optimality can be described by a single parameter 
Q > 0. This is easy to see by transforming to a coor- 
dinate system in which the chosen Gaussian is circular: 
in this coordinate system, Q is a function of the angle 
between the stable and unstable directions {Q ~ for 
a right angle and increases as the angle decreases to 0). 
In the coordinate system in which the manifolds are or- 
thogonal, Q will be a function of the aspect ratio of the 
Gaussian and its tilt relative to the axes. The general 
form of the short-time autocorrelation function is then 



^iin(t) = {a\a{t)) = 



\/cosh + iQ sinh At 



(25) 



As long as Q is not large, the qualitative scarring behav- 
ior is unchanged, although spectral oscillations become 
less pronounced and there is less scarring and antiscar- 
ring as Q increases. Analytic results are less easy to ob- 
tain for non-zero Q, but quantitative predictions can be 
readily produced for comparison with any experimental 
or numerical data. See lllitl for a discussion of the effect 



that non-zero Q has on antiscar strength, and thus on 
the probability to remain in a weakly open system. 

Finally, we conclude this section by mentioning that 
the analysis can be easily extended to test-states cen- 
tered off of the periodic orbit. In the coordinate system 
defined by Eq. |g, we may take a Gaussian centered at 
(qojPo), with width aVh in position space and apVK in 
momentum (where ap = l/cr): 



1 



1/4 



'(g— go) /2a- h+ipo{q—qo) 



(26) 



Then the autocorrelation function is given by pE 



exp 



sinh^ Xt/2 
cos 



X e 



V cosh Xt 

-i(pt — i{tEinh \t)qQPQ /h 



(27) 



We easily see that the linear recurrences, and thus scar- 
ring, are strong in a region of size crVh x apVh = h in 
phase space. Eq. ^ will be made use of in Section ^, 
where the tail of the wavefunction intensity distribution 
is obtained, summing contributions from test states lo- 
cated everywhere in phase space. 



IV. QUANTITATIVE TESTS 

We now provide numerica l ev idence supporting some 
of the conclusions of Section pill The data is taken from 



where the reader will find additional numerical tests 
along with a discussion of the generalized baker maps, 
the test systems used to obtain numerical data in this 
and the following section. The baker maps are a class 
of Bernoulli systems and a paradigm of hard chaotic be- 
havior. A symbolic dynamics allows one to easily identify 
the locations of the periodic orbits and to find their expo- 
nents. Furthermore, the local dynamics of these systems 
around any periodic orbit is exceedingly simple, the sta- 
ble and unstable manifolds always being oriented along 
the p and q directions, in accordance with Eq. ^. Several 
system parameters can be easily varied to collect wave- 
function statistics while keeping fixed the linearized dy- 
namics near a given periodic orbit. Semiclassical spectra 
and eigenstates can be efficiently computed, allowing for 
comparison between these and the exact quantum sta- 
tionary properties ||4C|l . 

In Fig. [ip appears a typical local line spectrum S{E) 
(Eq. ^ for a wavepacket centered on a periodic orbit 
with exponent A = 0.79. A total of 223 states exist in 
the energy interval [— 7r,7r]; a portion of this spectrum 
containing the peak of the linearized envelope Sn^iE) 
(dotted line) is shown in the figure. Two intermediate 
envelopes, defined by cutting off the sum in Eq. |l^ at 
T = 30 (a time long compared to the orbit period but 
short compared to the Heisenberg time) also appear in 
the figure. The solid curve is obtained from the exact 
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quantum calculation until |T| < 30, while the dashed 
curve comes from the corresponding semiclassical calcu- 
lation. We see that the large wavefunction intensities do 
appear in the energy region predicted by the linear enve- 
lope, and further refinement of the spectrum is obtained 
from intermediate-time dynamics. 

In Fig. I^b we again zoom in on the energy region in 
which the linearized envelope (now shown as a solid line) 
is peaked. The full spectrum S{E) is still marked by the 
vertical solid lines, while a semiclassical line spectrum 
appears as dashed lines. Though the quantum and semi- 
classical spectra differ in their detailed properties, they 
are both observed to fluctuate around the same linear 
envelope Sun- In Fig. |l|c we separate out these random 
fluctuations in the quantum spectrum from the short- 
time constraint 5*11,1, by dividing through all of the line 
heights by S\in- We see that the rescaled line heights dis- 
play energy-independent oscillations, equally strong near 
the peak and valley of the envelope. 

In |l^] the rescaled line heights of Fig. |l|c are collected 
for an ensemble of systems, and compared to the Porter- 
Thomas prediction of Gaussian random fluctuations. Be- 
cause we will be focusing on wavefunction intensities in 
the following section, we omit this data here, but will 
point out that the IPR for the spectrum of Fig. |l| is 7.93, 
which agrees well with the scar theory prediction 7.94 for 
this value of A when the random fluctuation factor F — 2 
is included (see Eq. |2^) . The IPR of the semiclassically 
evaluated spectrum for this system is 8.46. 

In Fig. g we plot the numerically obtained values 
of the IPR for a large number of systems, with the 
test wavepacket placed on orbits of varying exponent A 
(0.28 < A < 1.94). The actual value of the IPR is plot- 
ted versus the linear IPR prediction, < Sf^^ >, i.e. the 
amount of fluctuation in the linear envelope. The latter 
is obtained as in Eq. ^ by summing linearized return 
probabilities for a wavepacket centered on an orbit of ex- 
ponent A. Three sets of data points are shown in the 
figure (see 0]), and a line of slope = 2 is plotted for 
comparison. The data is consistent with the IPR being 
given by a product of the fluctuations in the linearized en- 
velope itself with Porter-Thomas fluctuations under the 
envelope. The spread of data points in the figure is also 
consistent with 0(l/\/^cff) fluctuations in the IPR from 
system to system, where N^g is the effective number of 
states under the linearized envelope [ pl| . 

In Fig. H we plot the 2-point correlation function 
< S{E + AE)S{E) > / < S{E) >2 of the spectrum 
shown in Fig. |l|a (plusses). In the same figure is plot- 
ted (with diamonds) the 2-point correlation function of 
the rescaled spectrum of Fig. ||c. Again, we see that all 
of the non-random structure in the spectrum arises from 
the presence of the linearized envelope Sn^, and that the 
rescaled spectral lines are completely uncorrelated. 



V. WAVEFUNCTION INTENSITY STATISTICS 

As an example of the calculations which can be eas- 
ily performed within the framework of scar theory, we 
discuss the tail of the wavefunction intensity distribu- 
tion, and its deviations from the exponential tail of RMT 
]l2| . This quantity has been analyzed previously for dis- 
ordered systems, where a log-normal tail is predicted l4l| ] 
and observed |Q. That analysis applies to scatterers 
that are small compared to the Fermi wavelength of the 
system. In the semiclassical limit, where the scatterer 
size is kept fixed while the wavelength is taken to zero, 
one obtains very different behavior, as we shall see be- 
low. This suggests that although RMT is known to be 
a good zeroth-order approximation for both quantum 
chaotic and diffusive behavior, the corrections to RMT 
may be qualitatively very different in the two kinds of 
systems. 

As we saw above in Section III, individual spectral 



lines in the local density of states S{E) show, in the 
semiclassical limit, chi-squared (Porter-Thomas) fluctu- 
ations around the linear envelope Sun- For complex 
eigenstates (in the absence of time-reversal symmetry) 
the chi-squared distribution has two degrees of freedom, 
and without scarring effects the probability of having a 
spectral line height greater than x is given by P{x) — 
exp(— a;). Here the intensity x is normalized to have a 
mean value of unity, i.e. Xn = iV|(n|a*'P") p, where 
\a'^°'P°) is the Gaussian wavepacket and N is the total 
number of states {N — 1/h for a classical phase space 
area of unity). Including scarring effects for an orbit of 
instability exponent A, this distribution is modified to 

P{qQ,Po,a,X,E,x) = exp{~x/Siiniqo,po,(J,X,E)) (28) 

for a eigenstate with energy E overlapping with a Gaus- 
sian centered at {qo,po) with width a. Thus, the aver- 
age intensity becomes 5'iin(9o,POj A, £*), and the entire 
distribution is stretched out by this factor. Now we per- 
form an energy average, remembering that for a map the 
quasi-energy E is defined to lie between and 2tt only. 
We notice that the tail of the intensity distribution will 
be dominated by the peak of the spectral envelope at the 
EBK quantization energy E — (p. Without loss of gen- 
erality we may set = 0, and perform the energy inte- 
gral using the saddle point approximation around E — 0, 
where Sun is maximized. We then obtain 

Piqo,Po,(T,\,x) = ^ J dEP{qo,PQ,a,X,E,x) 

1 ex.p{-x/Suniq,p,cr,X,E ^ 0)) 



2tt 



Sii„(0)^ dE-^ 



, (29) 



valid in the limit of large x. For small A, the sum over 
time steps can be replaced by an integral, and we have 
at go = Po = 



Sun{E) 



'iEt 



dt- 



\J cosh \t 



(30) 
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Now by dimensional analysis, 5*11,1(0) = Q/A and 
^Q^f^{0) = —W/X^, where Q and W are numerical con- 
stants. We thus obtain the h— independent tail of the 
intensity distribution for a wavepacket centered on a pe- 
riodic orbit, 



P{qo = 0,pa = 0, cr. A, x) = 



1 Q 



2tt VW 



XixX) 



(31) 



This is the total probability of having a wavefunction in- 
tensity enhanced by a factor x over the average intensity, 
on a periodic orbit of instability exponent A. The result 
is independent of cr, and thus of the aspect ratio of the 
test state in phase space. Notice that the exponential tail 
has been effectively stretched by a factor of Q/X, corre- 
sponding to the height of the peak of the linear envelope 
at small A. There is also a linear suppression factor of A, 
corresponding to the energy width of the peak in SuniE), 
and indicating that only a fraction scaling as A of all the 
eigenstates are effectively scarred. See the upper curve 
in Fig. ^ for a comparison of this result with numerical 
data. 

The region of validity of Eq. ^ is 



K A^i < a; < iV. 



(32) 



The first inequality ensures that many iterations of the 
periodic orbit contribute (so the sum over iterations can 
be replaced by an integral) and the scarring is strong. 
In fact, however, because of the large value of the nu- 
merical constant Q « 5.24, the formula works well even 
for exponents as large as log 2, as will be seen in the nu- 
merical study below. The second inequality says that we 
are in the tail of the distribution and the events are all 
coming from the peak of the linear envelope 5iin- The 
third inequality is a unitarity constraint: obviously our 
assumption of random fluctuations breaks down for in- 
tensities of order A'", when the entire wavefunction would 
be concentrated in a phase space area of order h. 

A similar analysis can be performed integrating over 
the phase space variables qo and po, to obtain the full dis- 
tribution of wavefunction intensities everywhere in phase 
space. The exponential exp (— a;/S'iin('i'07P07 • ■ ■)) must 
now be expanded to second order in qo,Po around the 
maximum qo — Po ^ 0- Then upon integration by the 
saddle point method we obtain a determinant factor fl^ ] 
of (ffy^-, again independent of cr. Z is another numer- 
ical constant, related to the falloff of Sun away from the 
fixed point qo = po = 0. The factor of h in front results 
from the fact that the dominant contribution to the tail 
comes from a region scaling as h surrounding the peri- 
odic orbit (compare Eq. |2^ and the discussion following) . 
Combining this determinant prefactor with the expres- 
sion in Eq. we obtain the tail of the distribution of 
overlap intensities after energy and phase space averag- 
ing. 



P{X,x) 



2 ZVW 



hX{xX)-^/^exp{-xX/Q). (33) 



This result is valid in the regime 

max(log TV, A"^) < a; < . 



(34) 



Here log N is the value of intesnity x near which the RMT 
exponential decay law exp(— x) reaches values of order 
h — 1/2ttN. Near this value of x, a crossover occurs be- 
tween the head of the distribution, which is dominated 
by non-scarred region of phase space and approaches the 
Porter-Thomas (RMT) prediction, and the tail, domi- 
nated by scarring, given by the expression above. The 
expression Eq. |3^ holds also for an ensemble of systems, 
all having one orbit with instability A. In principle, we 
should of course do a sum over all periodic orbits in 
a given system, however the tail will clearly always be 
dominated by the orbit with smallest A. Numerical data 
confirming Eq. 33 appears in the lower curve of Fig. ^. 

Finally, we consider an ensemble of systems where the 
value of the smallest exponent A varies from realization 
to realization, with distribution V{X) = CA" for small 
A. An example would be a billiard system with many 
movable disks, each of diameter much larger than the 
wavelength. Then using Eq. ^ and integrating over A 
we obtain 



n 



2 ZVW 



g"+i/2r(a + l/2)ax-(2+") . (35) 



Note that this is an uncontrolled approximation because 
we have integrated over A after having assumed xX to 
be large. However, if we had included higher-order cor- 
rections in (xX)^^ in Eq. ^ the scaling of P{x) would 
remain michanged, i.e. 



P{x) = Cf{a)hx~^^+"^ 



(36) 



with the dimensionless function f{a) somewhat differ- 
ent from that given in Eq. The important point is 
that the tail displays power-law behavior in the intensity 
X, a strong deviation from the exponential prediction of 
RMT. As with Eq. this asymptotic form is valid for 
values of x large compared to log TV and small compared 
to N. For small x we again expect a crossover to the 
Porter-Thomas form exp(— a;). For large x we expect a 
downward correction away from the a;^(^+") form, with a 
breakdown of the approximation occurring at some frac- 
tion of iV, depending on a. 

We now present numerical evidence for the results of 
this section, using the same ensemble of systems as in 
the numerical tests of the previous section. The data 
shown here was originally obtained in . In Fig. |^ we 
fix the instability exponent A of the least unstable peri- 
odic orbit at log 2, and obtain a cumulative wavefunction 
intensity distribution for a wavepacket centered on the 
orbit (upper thick curve), and a combined distribution 
for wavepackets randomly located over the entire phase 
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space (lower thick curve). The tails of the distributions 
compare weh with the scar theory predictions of Eq. |3l| 
(dashed curve) and Eq. |3^ (sohd curve). Notice that 
the generic wavepacket wavefunction intensity distribu- 
tion (lower thick curve) rolls over to the RMT prediction 
(dotted line) at small values of the intensity. The log- 
time in this system, near which value of x the crossover 
is expected to occur is Tiog ~ 7. 

In Fig. I the same distribution is shown after ensemble 
averaging over systems with classical orbits of different 
instability exponents. We see that at small values of x the 
data agrees well with the RMT prediction (dashed line), 
while the tail follows the power-law behavior of Eq. |35| 
(solid line). 



VI. IMPROVED (UNIVERSAL) MEASURES OF 
SCARRING 



We have noted previously an ambiguity in the choice 
of Gaussian wavepacket to serve as a test state for mea- 
suring scarring effects. Even after requiring that the 
wavepacket be optimally oriented relative to the invari- 
ant manifold directions, a one-parameter family of possi- 
ble Gaussians remains, parametrized by the aspect ratio 
(7^ (See Eq. ^. Together, these fill a hyperbohc region 
\qp\ < 0{h) encompassing the stable and unstable axes: 
the scarring phenomenon is present in this entire region 
and not just at the periodic orbit q — p — itself. [Com- 
pare this with the hyperbolic fringes in the Wigner spec- 
tral function obtained by Berry |^]. This fact suggests 
that the Husimi-space ambiguity may be resolved, and a 
better measure of scarring obtained, by including infor- 
mation about wavefunction behavior in this entire hyper- 
bolic region instead of choosing one Gaussian of arbitrary 
width a. 

A similar ambiguity arises when we consider a higher- 
period orbit of a map, or a periodic orbit of a continuous- 
time system. There, the test wavepacket is chosen to 
be centered at one of the many possible periodic points; 
again a less arbitrary measure of scarring should be ob- 
tainable by considering the wavefunction along the entire 
orbit. 

The obvious way to resolve these ambiguities is by av- 
eraging scarring intensity over wavepackets of different 
aspect ratios in phase space, or over wavepackets cen- 
tered at different periodic points, or both. Thus, we re- 
place the single wavepacket \a) by a density operator p, 
and define the scarring strength for a given eigenstate \n) 
as 



Inp = {n\p\n) 



(37) 



Then we can construct a wavepacket-averaged local den- 
sity of states analogous to Eq. |l|: 



5p(^) = ^/„p5(S-S„), 



(38) 



and a corresponding inverse participation ratio 

lPRp = Nj2lup- (39) 

n 

To obtain averaging over the wavepacket aspect ratio 
at a given periodic point, we write 



p^M / dte |a^gAt)(a„eAt| 



(40) 



The exponent A can be chosen to be the instability ex- 
ponent of the unstable orbit (Eq. but this choice is 
in fact arbitrary, and A can be reabsorbed into the over- 
all normalization factor J\f and the time cutoff Tq. In 
the absence of the cutoff Tq, the hyperbohc test state 
would be completely scale invariant, giving equal weight 
to Gaussians of all aspect ratios, from tall and thin to 
short and wide. The cutoff is, however, necessary be- 
cause the linearized dynamics of Eq. ^ is in fact only valid 
in a finite classical region around the periodic orbit (and 
also it eliminates possible normalization difficulties). To 
get an optimal measure of scarring, we choose the initial 
width a such that the wavepacket has the same aspect 
ratio in phase space as the linearizable region surround- 
ing the orbit in which the equations of motion are well 
described by Eq. ^. Then the initial wavepacket jao-) 
can expand by the same stretching factor along either of 
the two invariant directions before hitting the "edge" of 
the linearizable region. One easily sees that if A is the 
phase-space area of the linearizable region, we want to 
choose 



To 



c A 
— log — 
2A ^Ah 



(41) 



c is an arbitrary constant of order unity: it determines 
the exact suppression factor in Eq. Eo for a Gaussian 
just touching the boundary of the linearizable region. In 
the semiclassical limit, of course, A ^ h, and the pre- 
cise size A of the region needs to be determined only to 
within a multiplicative factor; it may indeed in this limit 
be taken as equal to the volume of the entire accessi- 
ble phase space. Notice that up to a constant, Tq has the 
same form as the log-time Tiog defined in Section [II . The 
Gaussian form of Eq. ^ in the time variable t is some- 
what arbitrary, and also not important in the h 
limit. See for a fuller discussion. 

To average over the position of the periodic point, we 
similarly may define a phase-space "tube" 



p^N dz\ 



(42) 



where the z coordinate parametrizes the periodic orbit in 
phase space, and at each periodic point the wavepacket is 
chosen to have width cr^ along the direction of the orbit 
and width a in the unstable direction at that point on 
the orbit. We recall that the Husimi function is nothing 
other than a Gaussian smoothing of the Wigner spectral 
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function; thus averaging the Husimi function around the 
periodic orbit coincides with the phase-space tube aver- 
aging of Fishman, Agam, et al. |^^. In |^^, the Gaussian 
width a is generally taken to be fixed as one goes around 
the tube; however that is a coordinate-dependent condi- 
tion, and in general one may take a to be an arbitrary 
function of z. 

More generally one can, of course, combine the averag- 
ing procedures of Eqs. to obtain a fully-universal 
measure of scarring. In the strong-scarring limit A <C 1 
one may also without loss of generality begin with one 
wavepacket only and allow the dynamics itself to gener- 
ate the density matrix p: 



Pdyn 



+00 

E 



|aiin(0)(aiin(i)l 



(43) 



Here |aiin(t)) is the original Gaussian \a) evolved in ac- 
cordance with the linearized dynamics: for example, if 
t is an integer multiple of the period Tp, then \a\in{t)) 
is centered at the same periodic point as |a), but with 
width (T cre-^*/^^. To is defined as before (Eq. |4l]), 
using the full instability exponent A for one iteration of 
the entire primitive orbit. X/P is the exponent per time 
step; hence the factor of Tp in Eq. 

To understand what gains may be achieved by the av- 
eraging process described here, we need to consider cor- 
relations in the local density of states between different 
wavepackets placed on the same periodic orbit. Indeed, 
these correlations will also be very important when we 
disc uss c onductances and phase space transport in Sec- 
tion VII. To begin, consider two wavepackets \a) and \b) 
and their overlap intensities Ina, Inb with the eigenstates 
I a) of the system. We may define a long-time averaged 
transport probability Pab p] as 



T-l 



P,,= hm i^|(a|C/*|6)| 

1 — !-00 1 



(44) 



For a nondegenerate spectrum we easily see 

Pab = |(a|n)|2 |(6|n)|2 = (45) 

n n 

In particular, the IPR corresponds to the special case 

IPRa - NPaa ^ N ^na ■ (46) 
n 

The Pab can be thought of as the covariance matrix of the 
local densities of states for different wavepackets, with 
Paa being the variances or diagonal matrix elements; the 
correlation between two local densities of states is then 
given by 



Cab 



V PaaPl 



(47) 



Consider first the case where \a) and \b) are exact time 
iterates of one another: \b) = |a(t)) for some time i. Then 
the two local densities of states are identical. 



Pab = Paa = Pbb - ^g{\) 



(48) 



and the correlation is unity. g(X) « tt/A for small A is the 
scarring IPR enhancement factor (Eq. |2^). On the other 
hand, if we consider two wavepackets that have the same 
short-time behavior but are otherwise unrelated (e.g. if 
I a) and |6) are located on different periodic orbits of the 
same length and instability exponent A), then the two 
have the same linearized spectral envelope but indepen- 
dent fluctuations under that envelope. Then 

Paa = Pbb = ^.g(A) 
Pab = ^5(A) 



Cab - - . 



(49) 



bb 



The correlation in this case is of order unity but still less 
than one. 

In the intermediate case, we may have two wavepackets 
one of which partly overlaps the other under short time 
evolution. Then the random spectral fluctuations under 
the two identical envelopes are partially correlated, and 
we get a correlation 1/F < Cab < 1 [Q- However, for 
small A, any two wavepackets optimally oriented along 
the same periodic orbit (whether centered at the same 
periodic point or not) will overlap each other with a very 
large amplitude under time evolution. The reason is that 
the wavepacket stretches very little each orbit iteration, 
and there will always be some n for which cre"^ is very 
close to any chosen width a' . Expressions for determin- 
ing the LDOS correlation Cab given A and the relative 
widths of the two wavepackets are given in [|l3[|_,_ where 
the predictions are also numerically tested [Fig. |6|. It is 
found that a scarring enhancement factor g{X) of 2 (cor- 
responding to instability exponent A ~ log 5), is already 
associated with a minimum correlation of 0.99 between 
the least correlated wavepackets on such an orbit. Or- 
bits that have smaller A and are thus more scarred have 
correlations Cab even closer to unity. 

In Fig. 1^ we plot the predicted and numerically ob- 
served correlations Cab between local densities of states 
for two wavepackets centered on the same periodic or- 
bit. The parameter z is defined such that wavepacket 
|a), under linear evolution, evolves to be centered at the 
same point as \b), but with width e'*'*^""'"^^ times that of 
|5), where n is an arbitrary integer. For integer z the two 
local densities of states are the same, and the correlation 
is unity. The smallest correlation is obtained for half- 
integer z. The upper and lower curves in the figure cor- 
respond to instability exponents A — log 5 and A — log 10, 
respectively. The data is obtained by averaging over an 
ensemble of kicked-baker maps |lJ , and the theoretical 
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curves are obtained from We notice that even for 

these very moderate instability exponents, the minimum 
value of Cab for the most uncorrelated wavepackets is al- 
ready very close to one. The correlations get even closer 
to one as A decreases. 

Thus in the strong scarring limit A ^ 1, there is in 
fact no ambiguity in choosing either the central point z 
or the aspect ratio of the test Gaussian on a periodic 
orbit. This also means that the density matrix averaging 
outlined above does not produce a scarring measure dif- 
ferent from that given by any one arbitrary wavepacket 
along the orbit. 

Yet it is possible to use our knowledge of the entire 
orbit and the linearized classical dynamics surrounding 
it to build an improved scarring test state. Beginning 
once again with a fixed point of a discrete-time map, we 
construct the coherent sum 



=M dte 



I' 



(50) 



(compare with Eq. |4C| ) . Again, Tq is a cutoff time (Eq. ^ ) 
for the linearized dynamics, and AA is a normalization 
constant. 9 is an arbitrary phase accumulated per time 
step. We notice that if the linear equations of motion 
(Eq. H) were exact, we could remove the cutoff Tq, and 
the state would be an exact eigenstate of the dynam- 
ics, with energy 4) + 9. The extent to which ^E" fails to be 
an exact eigenstate is determined by the nonlinear na- 
ture of the dynamics beyond the log-time To ~ | log?i|/A, 
whereas the failure of a single wavepacket \a) to be an 
eigenstate results already from the initial decay at the 
time scale 1/A. Just as scarring strengths and IPR's 
for a single wavepacket scale as 1/A because the auto- 
correlation function A^^ remains of order unity for the 
first 0(1/A) iterations of the orbit, so here the autocor- 
relation function should decay on a time scale of order 
To ~ I log?i|/A in the ?i — > limit, and scarring strengths 
should scale likewise. 

Indeed, in the strong-scarring semiclassical regime 
(A — > and ?i 0) we obtain |Q the short-time au- 
tocorrelation function for 



ALU) 



(51) 



Fourier transforming, we then have the spectral envelope 

St{E) = V2^roe-(^-^-^)'^«/2 , (52) 

centered at energy E = <j) + 9, with width scaling as 
A/I log?i|, and height and IPR [Fig.|^ scaling as | log?i|/A. 
In particular, the inverse particpiation ratio has the fol- 
lowing asymptotic form for small A and small h: 



IPR^, = F^Tq , 



(53) 



Note that the LDOS of a single wavepacket is always 
peaked at the "quantization energy" E = (f>, whereas 
by taking coherent superpositions we may construct an 



"off-resonance" state with a LDOS peaked away from 
this energy [see Fig. ||. For 9 — 0, the state \I' lives 
on the fixed point itself and on the invariant manifolds 
within the linearizable region around the orbit; for non- 
zero 9 it lives on hyperbolic regions asymptotic to the 
two manifolds [Fig. |^: 



qp 

h 



9/X 



(54) 



These hyperbolic structures living near the periodic or- 
bits but not on the orbits themselves have been stud- 
ied previously by Nonnenmacher and Voros [^ , and ob- 
served before that by Saraceno [^^ and by O'Connor and 
Tomsovic in the context of the baker map. The asso- 
ciation between energy shifts and hyperbolic structures 
surrounding the orbit is also evident in the phase of the 
dynamics in Eq. ^ for wavepackets centered off of the pe- 
riodic orbit. We see there that at times short compared 
to A^^, the phase accumulated per time step is shifted 
by the amount qoPoX/h (where {qo,Po) is the center of 
the wavepacket), and the peak of the spectral envelope 
in quasienergy space must therefore also be shifted by an 
amount of this order. 

The scarring measure given by Eq. ^ may truly be 
said to be universal and optimal, as it eliminates the am- 
biguity inherent in the Husimi space measure, and fully 
incorporates all knowledge not only about the location 
of the periodic orbit itself but also about the linearized 
dynamics around this orbit. We note here that the ex- 
pressions Eqs. are valid in the A — > limit, where 
scarring is strong. However, this improved measure of 
scarring can also be used for moderate A, where the ex- 
pressions (Eq. |5l| , |5^ ) need to be modified. Although the 
dependence of the spectral envelope on A is more compli- 
cated away from the A ^ 1 regime, the spectral envelope 
height and the IPR still scale as log?i in the semiclas- 
sical limit, for arbitrary A. Thus, even an unstable or- 
bit with moderate or large instability exponent, which 
will be rather weakly scarred using the usual Gaussian 
wavepacket measure, will always be strongly scarred us- 
ing the "universal scarmometer" ^, in the ?i — > limit. 
See [Q for a more thorough discussion of this point, as 
well as of finite Ti corrections to theoretical predictions of 
scar strength. 

We now briefly discuss the extension of the above 
methods to longer periodic orbits of maps and to contin- 
uous time. Here, even greater enhancement of scarring 
intensity is obtained because the universal test state ^' 
is naturally defined to live along the entire orbit. For a 
map of period Tp, we recall that there are Tp optimal 
energies 



E, = 



2TTk 



Tp 



. . . Tp - 1 



(55) 



that "quantize" the orbit in the EBK sense. The spec- 
tral envelope of a single Gaussian wavepacket will be 
peaked at all these Tp energies (see the discussion fol- 
lowing Eq. pj). We now choose one such quantization 
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energy and construct an optimal scarring test state ^ 
centered at energy Ek + 6/Tp: 



* =A/'^ J dte 



Here fp is a stretching factor, and (j)p is a phase, both 
defined by 



Xp/ 



(57) 



fp and (j)p take into account the fact that stretching and 
phase accumulation along the orbit may both be non- 
uniform; of course, /tp = 1 and (pTp — (j). The cor- 
rect phases 0p are particularly crucial for getting con- 
structive interference between different periodic points. 
Once again, 9 is an arbitrary off-resonance parameter. 
For 9 = 0, the test state lives on all Tp periodic points 
and their invariant manifolds; for 9 = 0(A), the test 
state occupies hyperbolic regions qp ^ h near these same 
periodic points. 

The short-time autocorrelation function and the spec- 
tral envelope take the same form as in the Tp = 1 fixed 
point case, making the substitutions 

To ^ ToTp , A ^ X/Tp , 9 ^ 9/Tp , (j) ^ Ek (58) 

throughout. Notice in particular that the relevant in- 
stability parameter is X/Tp, the exponent per time step, 
and not the exponent A per iteration of the entire or- 
bit which governs the strength of the single-wavepacket 
scarring measure. In particular, the peak height in the 
LDOS S^^{E) and the corresponding IPR both scale 
as (A/Tp)~^| log?i| for small X/Tp. Thus, as long as 
the Lyapunov exponent X/Tp remains small, the effects 
of longer orbits (up to the mixing time) can remain 
strong using the ^' test-state measure, whereas the single- 
wavepacket measure will fail to detect a strong signal 
from these longer orbits. For Tp > 1, a dramatic in- 
crease in scarring strength is observed even for moderate 
h, because of this additional Tp factor. 

For small A, the optimal test state can also be gen- 
erated dynamically: 



dyn 



t 



(59) 



(compare Eq. ^). Here \a) is a wavepacket of width a 
centered at any point along the periodic orbit. 

The expressions Eqs. also generalize in an ob- 

vious way to continuous time |p^ . Again, one chooses 
any energy E^ which quantizes the action in units of h, 
and an arbitrary energy offset h9/Tp if one wishes to 
measure off-resonance scars. The short-time autocorre- 
lation function again decays on the time scale TqTp ^ 
TpX-^ log A/h. In the energy domain it produces a sin- 
gle peak of height scaling as (Tp log^/?i)/A compared to 



the background spectral envelope (of width ?i/Tfi-eo, see 
discussion in Section HI) which enforces energy conserva- 
tion. This is to be compared with the ordinary scarring 
envelope which has peaks at all of the quantization en- 
ergies with height scaling simply as 1/A. Once again, 
we find that longer orbits can be detected as easily as 
short ones using the improved scarring measure^ and an 
additional logarithmic scarring enhancement factor is ob- 
tained in the semiclassical limit. 

We mention here the closely related work of de 
Polavieja, Borondo, and Benito on improved measures of 
scarring in Hamiltonian systems | |30| | and of Simonotti, 
Vergini, and Saraceno |Q on the surface of section in 
a billiard system. In both coherent superpo- 

sition of wavepackets is used to construct a better test 
state for measuring scarring effects. However, both these 
groups consider only one (arbitrarily shaped) wavepacket 
at each periodic point of the orbit; thus they obtain the 
enhancement factor Tp but not the extra |log?i| factor 
that comes from using the linearized dynamics to fill the 
entire hyperbolic region surrounding the orbit. See [T^ ] 
for a fuller discussion of this and other related issues. 

Numerical data in this section is obtained from an en- 
semble of kicked baker maps; see 13 1 for details. In 
Fig. 0, we show the IPR plotted as a function of the 
cutoff constant c (see Eq. ^), for various values of the 
system size N {N = 1/h is the number of Planck-sized 
cells in the classical phase space). Here the instability 
exponent for the period-one orbit is A = |log0.18|. The 
upper dashed curve is the N oo theoretical prediction, 
which converges to the asymptotic form of Eq. |5^ (lower 
dashed line) for large Tq (i.e. for exponentially small 
h, keeping fixed the cutoff parameter c). The linearized 
theory is expected to start breaking down for c > 1 (the 
rightmost six points on each data curve). 

In Fig. ^ are shown the smoothed (numerically ob- 
tained) local densities of states for off-resonance test 
states ^E*, for several off-resonance angles 9. The tallest 
peak corresponds to the on-resonance test state 6* 0; in 
the h, ^ limit the peaks would all be the same height 
and shape. The system size here is fixed at TV = 800, the 
cutoff parameter is c = 0.6, and once again A = | log 0.18| . 
In Fig. ^ are displayed Husimi plots of the universal test 
state * for (a) 9/X = 0.8 and (b) 9/X = 2.5. For = 
(not shown), 4* lives entirely on the orbit (located here in 
the center of the figure) and on the invariant manifolds 
(vertical and horizontal lines intersecting at the center 
of the figure). As 9 becomes non-zero, moves off of 
the manifolds and onto hyperbolic regions asymptotic to 
the manifolds (Eq. 54). Choosing the opposite sign for 



Of course, given a value of Ti only orbits of period less than 
the mixing time for that h can be profitably studied; our 
approach is based on the assumption Tp < Tiog, so that lin- 
earized dynamics remains valid over several iterations of the 
orbit. 
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9 would produce a test state occupying the other two 
quadrants. 

Finally, in Fig. |o| we show a numerically obtained, 
smoothed spectrum for a wavepacket placed on a period- 
2 orbit (double-peaked solid curve), and for the univer- 
sal test state ^ constructed on the same orbit (tall sin- 
gle peak). These agree well with the theoretical predic- 
tions (dashed curves) based on the linearized dynam- 
ics near the periodic orbit, which has a total exponent 
A = I log0.168| over the two-step period. One of two pos- 
sible on-resonance energies has been chosen for the test 
state '3/, which is again constructed using cutoff constant 
c = 0.6. The enhancement here is more dramatic than 
for the period-one orbits, due to the fact that universal 
scarring strength depends only on the exponent per unit 
time along the orbit, not on the orbit length itself. 



VII. OPEN SYSTEMS: SCARS AND ANTISCARS 

We now proceed to consider the last major topic to 
be covered in this review: the effect of unstable peri- 
odic orbits on weakly open systems, where experimen- 
tally accessible quantities such as resonance peaks and 
conductance curves can be measured. See for a more 
detailed exposition of the effect of periodic orbits on reso- 
nances and on the probability to remain in an open quan- 
tum chaotic system. A study of conductance curves in 
two- lead chaotic systems is now in progress |^ . 

Let -ffo be a quantum Hamiltonian for a classically 
chaotic (closed) system, and consider adding a small 
opening to the system, through which probability can 
escape. If the opening size corresponds to less than one 
open channel, so that quantum resonances will be non- 
overlapping, we can write an effective Hamiltonian for 
the open system: 



H = Ha-i^\a){a\. 



(60) 



Here |a) is a quantum channel associated with the open- 
ing, and F is the decay rate in that channel. \a) could 
be a Gaussian wavepacket at the location of the hole, or 
a position or momentum state. The main effect of the 
opening on a wavefunction \n) of the closed system is that 
it acquires a decay width proportional to the intensity of 
the wavefunction at the opening: 



r„ = r|(n|a)|2. (61) 

N\{n\a)\^ follow a chi- 
exp(— a;/2) for real 



In RMT, the intensities x„ 
squared distribution: P{x) = 

overlaps, and P{x) — exp (— a;) if the overlaps {n\a) are 
complex. On the other hand, if the decay channel \ a) hap- 
pens to be located on a short periodic orbit, we know that 
the distribution of decay rates is stretched by the factor 
S'g,-j(£^), this being the smoothed LDOS at \a). [Compare 
with the discussion of wavefunction intensity statistics in 



Section This will have a profound effect on the dis- 
tribution of resonance lifetimes in the system. Near the 
would-be EBK quantization energies where S'un > 1, res- 
onance lifetimes are decreased (by a factor scaling as 
for small A), while far from these energies resonances tend 
to be narrow compared with the expectations of RMT. 
In order to have a concrete quantity to focus on, let us 
consider the ensemble-averaged probability to remain in 
the system at long times. This quantity has been in- 
vestigated previously for disordered systems [^3| , and a 
log-normal long-time tail was found. We investigate here 
what deviations from RMT predictions can be obtained 
as a result of periodic orbit effects. 

First consider the RMT case, with complex eigen- 
states. Because mixing between the states of the closed 
system can be neglected in the small F regime, the total 
probability to remain in the system is given by a sum 
over these states: 



1 ^v-i ^oo 



i:Tt/N 



n=0 

oo 



1 



1 + Tt/N 



(62) 



Here N is the total number of states in the system; the 
classical decay rate is given by 



T/N 



(63) 



because only one channel has the possibility to decay. 
We see that at short times (t ^ F^j^), the probabihty to 
remain in the system is Prem(0 ~ 1 — Fdi, as expected, 
while at long times the power-law asymptotic behavior 



(64) 



is obtained. The analysis can be extended to M > 1 
open channels, where the power of the long-time behav- 
ior scales with M. Notice that at the classical level we 
expect an exponential decay law Prem(i) = exp(— Fdt), 
since the system is chaotic and the hole small enough that 
the probability distribution inside the system never de- 
viates from uniformity (the mixing time being very short 
compared to the decay time). 

Now take the open channel \a) to be centered on an 
unstable periodic orbit of exponent A. The intensities 
are then distributed in each energy range as a chi-squared 
distribution with mean Siin{E): 



P{x) 



1 



,-x/Sii„(B) 



(65) 



Then the probability to remain at long times is given by 

1 



1 1 



(66) 
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if initially only states with energy around E are popu- 
lated. Three qualitatively different energy regimes can 
be distinguished [Fig. At the quantization energy 

E = (j), S\in{E) reaches its peak; the peak height scales 
as A^^ for small A, and its width scales as A compared to 
the peak spacing. This is in fact the energy region which 
dominates the tail of the wavefunction intensity distribu- 
tion in Section^. Away from the peak, for |i? — 0| 3> A, 
the linearized LDOS falls off exponentially: 

S,,^{E) « . (67) 

A 

[Notice that for simplicity the above expression and those 
to follow are given in quasienergy units, where the spac- 
ing between successive scarring peaks is fixed at 27r. A 
factor of Tp/h, where Tp is the orbit period in real time 
units, needs to be inserted to convert E to a, real energy.] 
Within 0(A) of the antiscarring energy E — (j) + tt 
(this is the energy equidistant between two successive 
scarring peaks), the LDOS rolls over smoothly to attain 
the minimum value 

5n„(i? = + ^)«ye-'/'^ (68) 

The region within 0(A) of the antiscarring energy E — 
(j)+TT is thus responsible for producing the smallest wave- 
function intensities, and the narrowest resonances in the 
corresponding open system. This excess of exponentially 
small decay rates is as dramatic a signature of the under- 
lying classical behavior as the long wavefunction intensity 
tails found in Section ^ 

Energy-averaging allows us to sample resonances from 
all three energy regimes: wide resonances from near 
E — (j), very narrow resonances from near E — (p + i:, and 
intermediate- width resonances from the region described 
by Eq. |6^ (the same can be accomplished by varying a 
weak magnetic field at constant energy and thus sweep- 
ing through different values of the phase </>). The total 
probability to remain in the system at short times is given 
by 

Prcm = 1- < 5lin > Tt/N = 1 - T,it (69) 

(where we use the first line of Eq. ^ and notice that 

< •S'lin >= Aiin(O) = (a|a) = 1 by normalization). Thus 
at short times t <C T^^, the faster-decaying scarred states 
and slower-decaying antiscarred states always cancel ex- 
actly to produce a result consistent with our classical 
expectations. At long times, we use the second line of 
Eq. ^ to obtain 

Prc.„(t) = ^|4^. (70) 
t clt 

In the strong scarring regime of small A, the average 

< S^^ > is dominated by narrow resonances within 0(A) 
of the antiscarring energy E — cj) (Eq. |68| ), and gives an 
exponentially large (in A) enhancement over the predic- 
tions of RMT: 



Thus, the quantum mechanical probability to remain in 
the system at long times is very strongly dependent on 
the location of the opening relative to the unstable classi- 
cal orbits of the system. This is in contrast with the clas- 
sical long-time probability to remain, which for a small 
opening is independent of the location of the opening [the 
classical distribution continually being uniformly redis- 
tributed over the entire phase space, on a time scale very 
short compared to the decay time of the system]. See 
Figs. |l^,|l^ later in this section for a numerical example. 

Enhancements in the probability to remain at long 
times are obtained also for leads centered slightly away 
from the periodic orbit; in particular an enhancement of 
order \^e^ /^'^ is found as long as the lead is within a 
phase space area scaling as X^fi surrounding the orbit. 
Thus, the long-time probability to remain in the system 
averaged over all possible positions of the lead is given 

by 

[In principle, contributions from all the periodic orbits 
need to be added, however, if orbits with small A exist, 
they will clearly dominate any such sum.] The result is 
that at finite energy, exponentially large (in 1/A) devia- 
tions from RMT are found even after averaging over the 
lead position. In the ?i — > limit of any given classical 
system, the RMT behavior is recovered as the probabil- 
ity of the lead being within ?i of a short periodic orbit 
goes to zero. 

Classically, not only is the long-time probability to 
remain independent of the lead position, but it is also 
largely independent of the initial probability distribution 
in the system (as long as this distribution is not concen- 
trated in a very narrow corridor starting at the lead and 
having width scaling as the lead size). Similarly, the re- 
maining classical probability distribution at long times 
is always evenly distributed in the classical phase space, 
again with the exception of a similar narrow corridor. 
Both of these statements are modified by quantum me- 
chanical effects associated with the unstable periodic or- 
bits of the system |l^ . In the non-overlapping resonance 
regime of small F, the long-time probability to remain 
in state |6) given a decay channel at |a), as well as the 
total probability to remain in the system given an initial 
probability concentrated at \a) are both given by 

n 

If 1 6) is not located on the same short periodic orbit as 
I a) (and not related to it by symmetry), the intensities 
|(6|n)|^ follow a chi-squared distribution independent of 
the decay rates F„ ~ |(a|7i)|^, and the mean value of this 
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intensity at a given energy scales as S^^^, the smoothed 
LDOS at 1 6). Then at energy E we have 



SUE) 1 
SUE) Tcit 



(74) 



Again averaging over energy we obtain the long-time 
ratio of the remaining probability density at \b) com- 
pared to the average remaining probability density at 
long times: 



pb 



Pre 



(75) 



The ratio obtained above strongly deviates from unity 
if both I a) and \b) are located on short unstable periodic 
orbits, so that their linearized spectral envelopes S^^^ and 
"^Un nontrivial. Consider two specific cases. If S'j'Jj, 
and Sfijj are located on (distinct) periodic orbits of the 
same classical period, action, and stability exponent A, 
then the two envelopes are identical and we obtain the 
exponential suppression 



pb 

rem 
P 

^ rem 



=< [Sir' >-'^ 



(76) 



On the other hand, if the two spectral envelopes are ex- 
actly out of phase in the region of interest, i.e. S^^^{E) = 
Si^^{E + tt), the ratio is dominated by the peak in Sfj^, 
and becomes an enhancement factor A^^ for small A. We 
note that this suppression or enhancement of probabil- 
ity on orbits other than the one on which the lead is 
located is of course a quantum interference effect, having 
no counterpart in the classical dynamics of the system. 
Even stronger interference effects are observed on the or- 
bit where the lead is located, due to the very strong cor- 
relations (found in Section VI) between densities of states 
at different points on the same orbit; see |ll4[| f or details. 

The data in this section is taken from |]14|, and was 
originally obtained using an ensemble of perturbed cat 
maps. See JTif for a fuller discussion of the numerical 
evidence. In Fig. |l2| is shown the measured probabil- 
ity to remain in the system as a function of scaled time 
t' — Tc\t, for a generic lead location (plusses) and for a 
lead placed on an unstable periodic orbit with exponent 
A = 0.96 (squares). The classical exponential decay law 
(dotted curve) is shown for comparison. We see that the 
numerical results for a generically placed lead agree well 
with the RMT prediction of Eq. |6^, shown as a dashed 
line. The long-time probability to remain given a lead 
placed on the periodic orbit shows a large enhancement 
factor, consistent with the theoretical prediction of Eq. |70| 
(solid line). For this value of A, the predicted enhance- 
ment factor at long times is 11.0. The numerical data is 
obtained for a system size N — 120, and the decay rate 
per step in the exit channel is T — 0.1. 

In Fig. |l^, the predicted enhancement factor < S^^ > 
in the long-time probability to remain is plotted as a 
solid curve, as a function of the instability exponent A 



of the orbit on which the lead is located. The predic- 
tion is valid in the semiclassical limit iV — > oo; numerical 
data is presented for N = 120 (plusses) and N = 240 
(squares). We notice the very large enhancement fac- 
tors that are obtained even for moderate values of A, in 
accordance with the exponential asymptotic behavior of 
Eq. [T^. Finally, in Fig. ^ we give an example of the 
enhancement and suppression of intensity at long times 
found on an orbit other than the one on which the lead 
is located, as the relative phase between the two orbits is 
varied. In practice, a curve of this kind can be obtained 
by varying a magnetic flux line or a weak magnetic field 
enclosed by one of the orbits. The theoretical prediction 
(Eq. |7|) and numerical data are obtained for the case 
where both orbits have instability A = 1.76. Because the 
stretching factor for one iteration of the orbit is quite 
large, both the intensity suppression when the orbits are 
in phase, and the enhancement when the two orbits are 
out of phase are quite moderate in this example. 

An application of scarring methods to the study of 
conductance peaks in chaotic system with two leads is 
upcoming p5| ; here we summarize the ideas and the ma- 
jor results. We again consider the regime where the leads 
are narrow (and the temperature low), so that the reso- 
nance peak width is small compared to the level spacing. 
The height of the n-th conductance peak is then given in 
the appropriate units by 



G 



(77) 



where r„a , r„f, are the decay rates of the 7i-th resonance 
through each of the two leads. If the openings are of the 
same size, this reduces to 



G 



\{n\a)nn\b)\' 
|(n|a>P + |(n|6)|^ 



(78) 



where we have taken out an overall constant proportional 
to the size of the opening. In RMT, the two intensities 
|(n I a) p and |(n 1 6) p are of course taken to be independent 
chi-squared variables (of one or two degrees of freedom, 
for real or complex overlaps, respectively), producing the 
desired distribution of peak heights. Several situations 
are now considered which show marked deviations from 
the RMT predictions. 

(i) One of the leads, say a, is optimally placed on a 
short unstable periodic orbit, of exponent A. Then con- 
ductance peak heights show periodic fluctuations in en- 
ergy, with larger conductances at the scarring energies, 
and smaller conductances far from these energies. In par- 



ticular, for small A the mean value of 



becomes 



large (as A~^) at the peak scarring energies, and the peak 
heights are then governed entirely by the fluctuations in 
|(n|6)p. Near the minimum of S^^^, the effect is more 
dramatic: the peak height again becomes a simple chi- 
squared variable, here governed entirely by fluctuations 
in |(n|a)p, but the mean value of the peak height is expo- 
nentially small in A (see Eq. A fraction of 0(A) of all 
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peaks will have the suppression of Eq. this excess of 
exponentially small conductance peaks will of course per- 
sist after averaging over either energy or weak magnetic 
field strength. 

(ii) The two leads are placed on two orbits whose spec- 
tral envelopes are "in phase" in the energy range of in- 
terest. Say for simplicity that the two orbits also have 
the same exponent A, so S^^^ = S^^^. At each energy, 
the conductance peak heights now follow the same dis- 
tribution as in RMT, but the mean value is proportional 
to Sii^{E). Thus, at the scarring energies the conduc- 
tance peak heights are increased by 0(A~^) over the 
mean RMT value, whereas at the antiscarring energies 
the heights are suppressed by ©(A^e^'^ /^^). The energy- 
averaged (or magnetic-field averaged) distribution can 
also be obtained, and shows an enhancement in prob- 
ability for large and small peak heights (compare the tail 
of the wavefunction intensity distribution in Section ^ . 

(iii) The two orbits are "out of phase". Then strong 
suppression of conductance peak heights is predicted ev- 
erywhere, and the mean peak height is also exponentially 
small in A. The largest peak heights appear at energies 
intermediate between the maxima of Sy^^ and S^^^, where 

S'fijj. Even at these energies, the conductance is 



'-'lin 



exponentially suppressed compared to the RMT predic- 
tion. 

(iv) The two leads are located on the same orbit, or 
on two orbits related by a symmetry. Here |(n|a)p = 
and the peak height distribution is again given 
by a simple chi-squared distribution, with mean given by 
SiiniE). The energy-averaged mean peak height is inde- 
pendent of the stability exponent A, but once again we 
predict power-law (in A) enhancement of the conductance 
near the scarring energies, and exponential suppression 
near the antiscarring energies. 



the classical rate of escape away from these manifolds 
is slow. Examples include highly excited eigenstates in 
three-dimensional billiards and wavefunctions of few- 
body chaotic systems made up of identical particles with 
rotation and permutation symmetries |^6|. 

Vilela Mendes has analyzed "saddle scars," structures 
which arise from unstable harmonic motion along the sta- 
ble manifold of a classical saddle point Bluniel et 
al. extend quantum chaos techniques, including analy- 
sis of scars, to problems in which ray splitting surfaces 
are present |48|. Delande and Sornette study acoustic 
radiation from membranes, and find localization in the 
radiation directivity that can be attributed to wavefunc- 
tion scars 

The insights of scar theory can be extended also to 
short-time effects not associated with any classical struc- 
tures. These short-time recurrences may result from 
diffractive scattering or other hard quantum dynamics far 
from the semiclassical regime. Thus, in a diffusive sys- 
tem, the local, short-time behavior of a quantum particle 
may lead in certain models to strong effects on the statis- 
tics of anomalous wavefunctions, effects not captured by 
the standard field theory methods which generally 
ignore dynamics on scales shorter than a mean free path. 
The scar picture may prove to be a unifying framework 
for studying anomalous localization on graphs, lattices, 
in disordered systems, and in S-matrix eigenstates. These 
extensions of scar ideas to non-classical short-time behav- 
ior are now under development. 
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VIII. CONCLUSION 



We have seen that short unstable periodic orbits leave 
a strong imprint on the long-time and stationary proper- 
ties of a quantum chaotic system, even though the cor- 
responding classical dynamics loses all memory of these 
structures at long times. This insight can be extended 
to other situations where short-time classical structures 
lead to modifications of RMT predictions concerning 
the structure of quantum wavefunctions. Sundaram and 
Scharf study periodic orbits in complex phase space, 
with complex action, and show that these "ghost" orbits 
can also produce wavefunction scars. 

All of the discussion in the present paper applies to iso- 
lated orbits in a system of any dimension, though specific 
formulas need to be adjusted to take proper account of 
multi-dimensional Gaussian wavepackets and their short- 
time evolution. In higher-dimensional chaotic systems, 
classically invariant manifolds containing large families 
of individual orbits can also lead to scar- like effects, if 
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FIG. 1. In (a), the full spectrum is plotted along with 
the linear envelope (dotted line), an intermediate envelope 
corresponding to |r| < 30 (solid line), and a semiclassical in- 
termediate envelope (dashed line). In (b), a portion of the 
spectrum (solid) is compared to the spectrum obtained us- 
ing semiclassical eigenstates (dashed). In (c), the quantum 
spectrum has been divided out by the linear envelope of (a) 
(after 0). 



FIG. 2. A plot of the actual value of the inverse partic- 
ipation ratio (IPR) for a wavepacket centered on a periodic 
orbit (squares), versus the value predicted by the linear the- 
ory. IPR's in the semiclassical approximation are plotted us- 
ing the '-(-' symbol, and IPR's in a baker map with random 
matrix theory nonlinear behavior are plotted with triangles 
(after 0). 



FIG. 3. The two-point spectral correlation function of the 
scaled spectrum (diamonds) is compared to the correlation 
function of the raw spectrum (plusses), after ensemble and 
energy averaging (after [pT[). 
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FIG. 4. Cumulative wavefunction intensity distribution 

(a) as measured by a test state centered on a periodic orbit 
with instability A — log 2, plotted as the upper thick curve 
with scarring theory prediction given by dashed curve, and 

(b) averaged over the entire phase space of size 200/i, plotted 
as lower thick curve with theory given by solid curve. The 
dotted line is the Porter-Thomas law (after |[f^). 

FIG. 5. Cumulative wavefunction intensity distribution 
after ensemble averaging over systems with classical orbits of 
different instability exponents. The scar theory tail is given 
by the solid line, and the dotted curve is the RMT prediction 
(after jl|). 

FIG. 6. The correlation Cab (Eq. 0) between local densi- 
ties of states for two wavepackets |a) and \b) , lying on the same 
periodic orbit of instability exponent A, is plotted as a func- 
tion of width parameter z. The upper and lower curves cor- 
respond to A = log 5 and A = log 10, respectively. The width 
(along the unstable manifold) of wavepacket | b) is e^'"+^' 
times that of \a), where n is an integer (after |13|). 



FIG. 11. Smoothed local densities of states Siin(£') are 
plotted as a function of energy on a periodic orbit of instability 
exponent A = 0.20 (solid curve) and on an orbit with A = 0.15 
(dashed curve). We observe the peak at the EBK quantization 
energy _E = which scales as A~^, the exponential decay 
between £; = and = tt (Eq. ||), and the minimum at 
the anti-EBK energy E = n, which is exponentially small 
in A (Eq. ||). The RMT prediction S'iin(-B) = 1, which is 
applicable away from any short periodic orbit, is plotted as a 
dotted line (after |l4|). 

FIG. 12. The probability to remain in the open quan- 
tum system is plotted as a function of scaled time t' = Vdt. 
The classical prediction exp{—t') is shown as a dotted curve. 
The quantum probability to remain for a generic lead location 
(plusses) compares well with the RMT prediction 1/(1 + t') 
(dashed curve). For a lead placed on a short periodic orbit 
with instability exponent A — 0.96, we obtain the enhanced 
long-time probability to remain (squares), which agrees with 
the scar theory prediction < S^^ > /t' (solid line). The sys- 
tem size used for obtaining the data is = 120 and the decay 
rate per step in the exit channel is F = 0.1 (after fl^). 



FIG. 7. The inverse participation ratio (IPR) for hyper- 
bolic test state |\E') is plotted versus the log-time cutoff To (see 
Eqs. po||4l|), for various values of system size A^. From bottom 
to top, the five curves correspond to A'^ = 50, 100, 200, 400, 
and 800. For each A'^, 26 points are plotted, for c = (1.1)'', 
j = —20... + 5. The orbit has exponent A — |log0.18|. 
The upper dashed curve is the N ^ ex theoretical predic- 
tion, which converges to the asymptotic prediction of Eq. ^ 
(lower dashed line) for large Tq. The linearized theory is ex- 
pected to start breaking down for c > 1 (rightmost six points 
on each data curve) (after |^). 

FIG. 8. Smoothed local densities of states are plotted 
for the off-resonance universal hyperbolic test state l^t), for 
off-resonance angle 8 ranging from (tallest peak) in steps 
of 7r/10 (to the left), through —2n, on a periodic orbit with 
exponent A — \ logO.lSj. The system size is A'^ = 800, and the 
cut-off constant c is set to 0.6 (see Eq. |4l| ) (after p^). 

FIG. 9. Husimi plots of the universal test state ^' 
for off-resonance parameter values (a) 0/X = 0.8 and (b) 
9/X — 2.5. The linearizable region is taken to be much larger 
than the displayed area of size 12v^ x 12%/^, centered on the 
periodic orbit (after psf). 



FIG. 13. The long-time enhancement factor of the prob- 
ability to remain in a system when the lead is placed on a 
periodic orbit is plotted as a function of the instability ex- 
ponent of the orbit. Data is shown for A' — 120 (plusses) 
and TV — 240 (squares). The A'' — > oo theoretical prediction 
< S^^ > is shown as a solid curve. We see the exponen- 
tial increase in the probability to remain as the exponent A 
decreases (the A ^ asymptotic form is given in Eq. [7l[ ). 
For large A, the enhancement factor converges to 1, the RMT 
prediction (after |l4t ). 

FIG. 14. The predicted relative intensity at long times 
on a periodic orbit other than the one containing the lead is 
plotted as a function of the relative phase between the two 
orbits (solid curve). For reference, the phase-space averaged 
intensity is plotted as a dashed line. Intensity suppression 
is predicted and observed when the orbits are in phase, and 
enhancement is seen when the orbits are exactly out of phase. 
The data is collected for an ensemble of systems with A'^ — 80, 
and each orbit has instability exponent A = 1.76 (after |l^). 



FIG. 10. Smoothed local densities of states for a Gaussian 
wavepacket placed on a period-2 orbit (double-peaked solid 
curve), and for the universal test state 'I' constructed on the 
same orbit (tall single peak). The dashed curves represent 
theoretical predictions based on the linearized dynamics near 
the periodic orbit in question. The periodic orbit has a total 
exponent A = | log0.168| over the two-step period (after lla]). 
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